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ABSTRACT 


A  direction  finding  antenna  based  on  residue  number  systems  (RNSs)  is  presented. 
Special  spacing  requirements  for  the  array  elements  are  derived  and  processing  of  the  array 
ouq)ut  for  high  speed  direction  of  arrival  (DO  A)  estimates  is  discussed. 

The  RNS  antenna  processor  encodes  the  interferometer  phase  response,  which  is  a 
sawtooth  folding  waveform.  By  design,  the  phase  response  of  each  element  pair  folds  with 
folding  period  equal  to  the  chosen  modulus. 

If  DOA  samples  are  generated  by  stepping  the  emitter  direction  between  -90  and  90 
degree  in  small  increments,  some  samples  will  fall  about  the  code  transition  points  and 
result  in  encoding  errors.  The  resulting  DOA  estimates  firom  such  a  realizable  system  could 
contain  large  spikes  or  “glitches”  at  these  points.  These  encoding  errors  in  tire  resolved 
DOA  can  be  reduced  by  interpolation.  Three  primary  methods  will  be  discussed  to  compare 
the  capability  of  removing  the  glitches:  LSB-Shift  Method,  Random-LSB  Method,  and 
Shift  Last  Good-Sample  Method.  A  comparison  of  the  performance  of  the  three  methods  is 
made  on  the  basis  of  simulation  data. 
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INTRODUCTION 


A.  OVERVIEW 

Direction  finding  (DF)  is  the  area  of  electronic-warfare  support  (ES)  whose 
objective  is  to  find  the  direction  of  arrival  (DO A)  of  a  targeted  transmitter.  ESM  performs 
search,  interception,  location,  and  identification  of  hostile  radiations.  DF  systems  passively 
examine  the  spectrum  use  by  hostile  emitters  and  process  the  signals  to  obtain  enemy  DOA 
information. 

A  general  assumption  in  most  DF  systems  is  that  the  received  field  exhibits  far-field 
plane-wave  behavior  with  linear  polarization.  In  reality,  incident  fields  are  often  non-planar 
with  phase-front  distortion  caused  by  multipath  and  scattering.  Typical  DF  system 
techniques  assume  that  only  a  single  signal  source  is  received.  Because  of  multipath 
ionospheric  propagation  (skywave),  and  wave  scattering  by  objects  close  to  the  DF  site  or 
co-channel  transmitting  stations,  signals  via  multiple  paths  and  sources  are  often  received  at 
DF  sites.  A  DF  system  designed  to  handle  one  source  at  a  time  can  work  in  a  multiple 
source  environment  if  other  signals  are  well  separated  in  frequency,  time,  or  direction  of 
arrival. 

In  practice,  DF  stations  often  receive  signals  from  several  sources  in  the  same 
frequency  range  at  the  same  time,  and  from  many  different  angles  of  arrival  [Ref  1].  To 
resolve  multicomponent  wavefields,  modem  techniques  of  spectral  estimation  methods  are 
used.  These  methods  include  eigenanalysis  and  linear  prediction. 

Modem  direction  finding  systems  fall  into  three  main  categories:  interferometer- 
techniques,  scanning  beam,  and  simultaneous-multiple  beam.  Interferometer  systems  based 
on  phase-comparison  techniques,  have  the  advantage  of  a  fast  response,  but  generally  use 
wide  coverage  antennas  which  result  in  low  sensitivity.  In  addition,  they  require  relatively 
complex  microwave  circuitry  that  must  maintain  a  precise  phase  match  over  a  wide 
frequency  band  under  extreme  environmental  conditions.  When  high  accuracy  is  required, 
wide  baseline  interferometers  are  used  with  ambiguity  resolving  circuitry.  The  mechanically 
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scanning  beam  requires  only  a  single  antenna.  However,  this  techmque  has  the  disadvantage 
of  a  low  probability  of  intercept.  The  simultaneous  multiple-beam  system  uses  an  antenna 
or  several  antennas,  forming  a  number  of  simultaneous  beams,  thereby  retaining  the  high 
sensitivity  of  the  scanning  beam  approach  while  providing  fast  response.  However,  it 
requires  many  parallel  receiving  channels,  each  with  full  frequency  coverage. 

This  thesis  is  concerned  with  a  technique  for  radiating  element  placement  and  signal 
processing  based  on  the  residue  number  systems  (RNS).  In  principle,  the  direction  of  arrival 
can  vary  in  azimuth  (0-360  degrees)  and  elevation  (0-90  degrees).  However,  this  paper 
deals  exclusively  with  the  azimuth  angle  of  arrival  for  a  linear  array.  In  addition,  an 
underlying  assumption  is  that  the  receiver  is  far  enough  from  the  transmitter  so  that  the 
arriving  wavefront  is  essentially  planar. 

The  basic  concept  of  residue  arrays'  was  presented  in  a  previous  thesis  by  Rodrigues 
[Ref  2].  A  simple  block  diagram  of  such  a  system  is  shown  in  Figure  1.  It  is  essentially  an 
interferometer  type  DF  system  consisting  of  IV- 1  pairs  of  elements.  Each  pair  includes 
element  1  so  that  there  are  iV-1  pairs  if  the  total  number  of  elements  is  N.  The  number  of 
pairs  is  determined  by  the  number  of  moduli  employed  in  processing  the  signal.  It  will  be 
seen  that  the  number  of  moduli  controls  the  angle  resolution,  the  element  spacings,  and  the 
number  of  comparators  required  in  the  processor.  The  phase  difference  between  each  pair 
of  elements  is  obtained  using  a  phase  detector.  In  principle,  phase  or  amplitude  comparison 
can  be  used,  but  Rodrigues  has  shown  that  the  phase  is  less  sensitive  to  variations  in  the 
DOA  and  element  factors.  The  phase  differences  are  quantized  and  encoded  using  the  RNS 
[Ref  3].  The  encoded  data  from  all  pairs  is  combined  to  obtain  a  DOA  estimate. 

This  DF  method  has  several  attractive  features.  First,  arbitrarily  fine  unambiguous 
angle  resolution  can  be  achieved  with  as  few  as  three  elements.  Secondly,  the  method  is 
fast,  providing  a  DOA  estimate  limited  only  by  the  RNS  processor  speed.  The  processor 


‘  This  term  will  be  used  throughout  the  thesis  to  refer  to  arrays  with  element  locations  or 
encoding  algorithms  based  on  residue  or  symmetrical  number  systems. 
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employed  here  is  based  on  parallel  digitization  of  the  N-l  phase  signals,  and  therefore 
nearly  instantaneous  DOA  estimates  are  obtained. 


DOA 


Figure  1.  DF  Array  with  RNS  Processing 

The  RNS  DF  method  has  several  potential  problems.  They  are  primarily 
performance  limitations  imposed  by  imperfections  in  the  system  hardware.  Several  of  these 
occur  in  the  array  beamforming  network: 

1 .  Mutual  coupling  variations  due  to  the  non-periodic  separation  of  elements. 

2.  Imbalance  in  the  weights  and  transmission  lines. 
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3.  Low  accuracy  of  the  phase  comparators  due  to  low  signal  levels. 

Other  errors  are  introduced  in  the  RNS  processing: 

1 .  Encoding  errors  that  arise  from  variations  in  the  comparator  thresholds. 

2.  Increase  in  quantization  error  due  to  frequency  deviations. 

The  most  serious  problem  is  the  encoding  errors  which,  because  of  the  non-sequential 
binary  encoder  outputs,  can  cause  extremely  large  errors  in  the  DOA  estimate.  In  fact,  the 
net  effect  of  all  of  the  errors  listed  above  can  be  reduced  to  an  equivalent  encoding  error. 

This  thesis  examines  three  interpolation  schemes  that  can  be  used  to  reduce  the 
encoding  errors  which  result  from  quantization  in  the  RNS  signal  processing.  Simulation 
data  shows,  for  example,  that  small  errors  in  the  comparator  thresholds  can  result  in  large 
DOA  errors.  Interpolation  methods  can  be  used  to  reduce  the  occurrence  and  magnitude  of 
the  errors. 

Three  types  of  interpolation  are  investigated:  (1)  a  shifted  least  significant  bit  (LSB) 
method,  (2)  a  randomized  LSB  method,  and  (3)  a  shift  last  good  sample  method.  The 
effectiveness  of  these  three  approaches  is  examined  using  computer  simulation. 

B.  THESIS  OUTLINE 

Chapter  II  summarizes  common  direction  finding  techniques  such  as  amplitude 
response,  phase  difference,  and  time  delay  methods.  Chapter  III  covers  array  beamforming 
and  the  mathematical  development  for  a  linear  array  anterma.  Also,  direction  finding  for 
multiple  emitter  DOAs  is  introduced  by  using  the  array  correlation  matrix.  Chapter  IV  is 
devoted  to  the  process  of  encoding  the  phase  response  with  the  RNS.  First  the  RNS 
encoding  procedure  is  briefly  introduced.  The  RNS  antenna  architecture  is  then  derived 
with  simulation  results  shown  to  demonstrate  the  array  performance.  Next,  we  examine 
encoding  the  phase  response  with  the  RNS.  The  RNS  encoding  procedure  is  briefly 
described  and  simulation  results  are  presented.  In  Chapter  V,  the  encoding  error  correction 
algorithms  are  derived.  Three  primary  methods  are  investigated  to  remove  the  quantizing 
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errors  or  “glitches.”  A  comparison  of  the  performance  of  the  three  methods  is  made  on  the 
basis  of  simulation  data.  Chapter  VI  presents  a  summary,  some  concluding  remarks,  and 
recommendations  for  future  work. 
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n.  DIRECTION  FINDING  PRINCIPLES  AND  METHODS 


As  early  as  1899  it  was  recognized  that  electromagnetic  waves  could  be  aimed  or 
beamed  in  a  specific  direction  by  the  use  of  antennas  such  as  loops  or  spaced  dipoles.  For 
receiving  purposes  this  meant  that  the  direction  of  arrival  of  a  radio  wave  could  be 
determined  by  moving  or  steering  a  directive  anteima  beam  in  azimuth  until  a  maximum 
signal  was  obtained.  The  beam  could  be  steered  by  physically  moving  an  antenna  or  by 
appropriately  summing  the  outputs  of  the  individual  elements  of  an  array. 

The  science  of  direction  finding  has  been  the  subject  of  much  study  since  the  early 
1900’s.  A  major  improvement  over  early  DF  systems  came  wifii  the  development  of  the 
Adcock  DF  system  employing  four  or  more  spaced  elements  with  an  overall  aperture  of 
one-half  wavelength  or  less.  The  Adcock  system  achieves  an  instrumental  accuracy  of  about 
5  degrees. 

Since  the  end  of  World  War  D,  there  has  been  a  trend  toward  larger  aperture 
systems.  For  example,  Doppler  DF  systems,  in  which  fi:equency  modulation  (FM)  by  a 
sinusoid  is  induced  on  the  received  signal  by  a  rotating  antenna  or  its  equivalent,  usually 
range  in  size  from  one-half  to  five  wavelengths  and  are  capable  of  accuracies  of  2  degrees 
or  less  [Ref.  4]. 

The  Wullenweber  DF  system,  which  was  developed  in  Germany  in  the  early  1940's, 
consists  of  a  large  circular  array  on  the  order  of  ten  wavelengths  in  diameter  [Ref  5].  At  any 
instant  of  time  only  one  segment  of  the  array  is  used,  and  with  proper  phasing  of  the 
elements  in  use,  this  segment  is  made  equivalent  to  a  planar  array.  By  commutation  the 
planar  array  is  in  essence  caused  to  move  in  a  circle.  An  accuracy  of  0.5  degree  or  less  is 
possible  with  this  system. 

At  microwave  frequencies,  antennas  having  aperture  sizes  of  several  wavelengths 
and  with  extremely  directive  patterns  can  be  physically  rotated  to  give  accurate  DOA 
information.  An  alternate  method  which  is  very  accurate  for  large  values  of  signal-to-noise 
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ratio  (SNR)  is  the  amplitude-sensing  monopulse  technique  in  which  two  identical  directive 
antennas  are  placed  some  distance  apart. 

A.  SPATIAL  DIRECTION-FINDING  PRINCIPLES 

The  main  function  of  a  direction  finder  is  to  find  out  the  direction  of  arrival  of  an 
incident  electromagnetic  wave  relative  to  the  coordinates  of  the  direction  finding  site. 
Figure  2  shows  a  representative  direction-finding  spatial  coordinate  system  with  the  DF 
located  at  the  origin  coordinate  point.  The  direction  of  arrival  is  specified  by  the  azimuth 
angle,  measured  from  the  y-axis,  and  the  elevation  angle  measured  fi-om  the  z-axis. 

General  direction-finding  principles  assume  that  the  electromagnetic  field  incident 
on  the  DF  exhibits  a  far-field  planar  wave  structure  with  linear  polarization.  As  depicted  in 
Figure  2,  the  radial  E-field  is  zero.  The  E-field,  and  the  H-field  are  in  space  quadrature 
[Ref  6]. 


The  direction  of  propagation  is  indicated  by  Poynting's  vector,  P .  In  practice,  the  incident 
electromagnetic  field  is  usually  nonplanar  with  phase-front  distortion,  created  by  multipath 
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and  scattering,  and  depolarization  produced  by  nonuniform  ionospheric  propagation 
effects.  A  full-capacity  generic  DF  should  measure  direction  of  arrival  in  three-dimensional 
space.  However,  many  operational  situations  require  only  measurement  of  the  directional 
component  in  the  azimuth  plane. 

The  measured  angle  of  arrival,  is  usually  called  a  bearing  angle  or  line  of 
bearing  (LOB).  Before  examining  residue  arrays,  the  capabilities  of  several  conventional 
DF  methods  are  presented. 

B.  DIRECTION  FINDING  TECHNIQUES 

1.  Overview 

Common  direction-finders  determine  direction-of-arrival  information  by  one  of 
three  measurement  methods  [Ref.  6]: 

1.  Amplitude  response, 

2.  Phase  difference, 

3.  Time  delay. 

The  angle  of  arrival  is  converted  into  a  voltage  analog  using  amplitude  response,  phase 
difference,  or  time  delay.  Basic  conversion  techniques  are  presented  in  the  following 
sections. 

2.  Amplitude  Response 

Two  t5q)es  of  amplitude  response  methods  are  used  to  obtain  DF  mformation:  direct 
and  comparative.  For  direct  amplitude  response,  the  directional  properties  of  elemental  H- 
field  loops  and  E-field  dipoles  provide  response  minima  as  they  are  rotated  in  azimuth. 
Bearing  data  are  acquired  at  these  minima  response  orientations. 

Elemental  loops  and  dipoles  exhibit  figure-eight  responses  as  depicted  in  Figure  3 
for  the  azimuth  plane.  Response  maxima  are  broad  but  the  minima  (nulls)  are  sharp.  The 
amplitude  response  derivative  is  greatest  around  the  null  response  position.  Therefore,  a 
dynamic  measure  of  null  position  provides  the  best  bearing  estimates. 
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Sequential  amplitude  measurements  are  made  as  the  directive  pattern  is  rotated  in 
aTi’rrmth  and  the  measured  bearing  is  indexed  to  the  azimuth  at  which  some  predetermined 
amplitude-null  property  occurs. 


Figure  3.  Azimuth  Plane  Response  Patterns 

The  horizontal  E-field  dipole  is  rarely  used  for  DF  purposes,  primarily  because  most 
signals  of  interest  are  vertically  polarized  and  of  orthogonal  polarization  to  the  horizontal 
dipole  [Ref.  6]. 

For  comparative  amplitude  DF,  multiple  antennas  are  used  to  obtain  orthogonal 
patterns.  Bearing  information  is  obtained  fi:om  the  ratio  of  the  amplitudes  of  the  orthogonal 
patterns.  Comparative  amplitude  DF  techniques  provide  instantaneous  bearing  information. 

3.  Phase  Difference 

Bearing  measurements  using  a  phase  difference  require  at  least  two  separated 
antennas.  A  plane  wave,  arriving  at  an  angle  other  than  the  normal  to  the  baseline  between 
the  two  elements,  arrives  at  one  element  earlier  than  the  other.  The  time  lag  between  the 
antennas  produces  a  RF  phase  delay  or  a  differential  RF  phase  between  the  antenna  outputs. 

Figure  4  illustrates  the  basic  phase-delay  technique.  An  incident  plane  wave  arrives 

at  an  incident  angle  9  at  antenna  1  inducing  a  voltage  exp(70r) .  After  propagating  the 
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distance  d  sinO ,  the  incident  plane  wave  induces  a  voltage  exp(7<» t-  0^  in  antenna  2, 
where  the  phase  delay  given  by 

0 = (Ind  /  X)  sin0  .  (1) 


Figure  4.  Phase-difference  DF  Parameters 

Therefore  the  bearing  angle  is  encoded  as  a  function  of  phase  delay.  For  the  phase 
difference  technique,  the  bearing  angle  is  computed  by  using  Eq.  1,  where  phase  difference 
is  measured  and  baseline  distance  d  and  wavelength  2.  are  known.  This  phase-delay 
technique  can  experience  phase  and  bearing  ambiguities.  If  d>X/2  phase  ambiguities 
exist;  however,  bearing  ambiguities  are  resolved  by  using  auxiliary  antennas  to  establish 
known  phase  references  [Ref  6]. 

A  variant  of  the  phase  measurement  method  is  Doppler  direction  finding.  The 
received  frequency  fi:om  a  moving  antenna  experiences  a  Doppler  shift.  If  the  antenna 
moves  along  a  radial  to  the  emitter  position,  the  Doppler  shift  is  maximized.  No  Doppler 
shift  occurs  if  the  antenna  movement  is  tangential  to  the  direction  of  propagation.  These 
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Tnavima  and  minima  Doppler  shifts  may  be  used  to  measure  die  bearing  of  the  incident 


signal. 

4.  Time  Delay 

An  electromagnetic  wave  propagating  between  two  separated  antennas  experiences 
a  time  delay.  This  time  differential  of  arrival  (TDOA)  contains  bearing  angle  information. 
In  Figure  5,  the  differential  time  delay,  At ,  between  antennas  1  and  2  is  given  by 


(2) 


where: 

c?  =  the  distance  between  antennas  1  and  2  (meters); 
c  =  velocity  of  light; 

0  =  bearing  angle  (degrees). 

Historically,  TDOA  has  been  used  for  pulsed  emissions  such  as  radar.  Also,  relatively  long, 
multiple-wavelength  baselines  are  commonly  used  [Ref  6]. 
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in.  ARRAY  BEAMFORMING  THEORY 


A.  OVERVIEW 

Direction  finding  systems  consists  of  an  antenna  (or  collection  of  antennas), 
beamforming  networks,  detectors,  and  signal  processing  equipment  It  operates  by  receiving 
and  then  analyzing  the  nature  of  the  emitter  signal.  In  noise-limited  situations,  the  source 
signals  may  consist  of  a  self-product  and  additive  receiver  noise.  In  more  complicated 
situations,  it  may  also  contain  a  noise  component  due  to  backscatter  from  the  surrounding 
environment  (clutter).  A  further  soxirce  of  complication  may  be  the  presence  of  interferences 
produced  by  jammers  beamed  at  the  DF  site  from  directions  unknown  to  operators  at  the 
site.  In  the  following  analysis,  unless  otherwise  noted,  these  complications  are  neglected. 
Emitter  signals  are  modeled  as  time-harmonic  plane  waves,  narrowband  assumptions  are 
made,  and  noise  is  assumed  to  be  white  Gaussian.  Before  examining  residue  arrays,  the 
capabilities  of  several  conventional  DF  methods  are  presented. 

B.  BASIC  SIGNAL  AND  NOISE  MODELS 

Figure  6  shows  a  uniformly  spaced  linear  array  consisting  of  N  elements,  and  with 
M  plane  waves  arriving  at  the  array  from  distinct  directions.  We  assume  that  all  these  plane 
waves  are  narrowband  with  the  same  carrier  frequency^  •  The  noiseless  signal  produced  at 
the  nth  element  of  the  array  by  the  wave  from  source  m,  denoted  as  the  wth  plane  wave,  is 
as  follows: 


s{n,m,t)  =  cos 


2nfct  + 


N  +  l 


'  A  I  7 


where: 


(3) 
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n  =  1,2,. ..,N  ( element  index) 
m  =  1,2. . . ,M  ( emitter  index) 

=  carrier  frequency 

-  .  .  2;r  j- 

=  vector  wavenumber  of  the  Mh  incident  plane  wave  =  — 

/L 

z  =  unit  vector  along  the  line  of  the  array 
t  =time 

The  parameter  Am  denotes  the  amplitude  of  the  signal  s[n,m,t). 


1  2  3  •  •  •  N-1  N. 

Figure  6.  Geometry  of  an  iV  Element  Linear  Array  with  M  Incident  Plane  Waves 

Let  Gm  denote  the  angle  of  arrival  of  the  mth  plane  wave,  measured  with  respect  to 
the  normal  to  the  array,  as  illustrated  in  Figure  6.  Then  we  may  express  the  dot  product  of 

the  two  vectors  and  z  as  follows  [Ref.  7]: 
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(4) 


^  liid  .  - 

where  d  is  the  interelement  spacing  and  X  is  the  received  radar  wavelength.  Accordingly,  we 
may  rewrite  Eq.  (3)  as: 

iv+n 


s{n,m,i)=  ^cos 


2Tt/j  +  2n^(^n 


2  J 


sin  6. 


(5) 


We  may  simplify  the  representation  of  the  signal  s{n,m,t)  by  adopting  the  complex 
notation  for  narrowband  signals.  In  particular,  we  may  represent  the  signal  s{n,m,t)  by  its 
complex  amplitude,  defined  by: 


5  («,/«)  =  4,  exp 


N  +  \' 


n- 


V 


sin0. 


(6) 


Now  define  the  electrical  phase  angle  from  element  to  element  along  the  array  as: 


d  .  ^ 

<t>m  =  27i-sme^  • 


(7) 


Also,  let  am  denote  the  complex  amplitude  of  the  signal  s[n,m,t)  measured  at  the  center  of 
the  array;  that  is: 


(N  +  \ 


-,m 


(8) 


n=- 


Note  that  if  the  number  of  elements  N  is  odd,  the  center  of  the  array  represented  by 
^  corresponds  to  the  center  element.  On  the  other  hand,  when  N  is  even,  the  center 


of  the  array  is  a  fictitious  point.  In  any  event,  we  may  use  Eq.  7  to  rewrite  Eq.  6  in  the 
simplified  form 


5(«,m)  =  a„exp 


j\n 


N  +  \ 
2  J 


(9) 


The  signal  measured  at  the  output  of  any  element  in  the  array  differs  from  the  signal 
actually  received  by  that  element  by  an  amount  attributed  to  noise.  In  practice,  this  noise  is 
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modeled  as  a  white  ergodic  random  process.  To  limit  the  effects  of  the  noise,  it  is  customary 
to  filter  the  ouQ)ut  of  each  element  in  the  array  to  the  same  narrow  band  of  fi-equencies 
occupied  by  the  actual  received  signal.  As  a  result  of  this  operation,  we  have  an  observed 
gignal^  which  for  each  element  of  the  array,  consists  of  the  actual  received  signal  plus  a 
narrowband  noise.  For  element  n  we  may  describe  this  narrowband  noise  by 

w(n,0=£,cos(2n//  +  |5„)  (10) 

where  the  amplitude  Bn  is  a  Rayleigh  distributed  random  variable  and  the  phase  is 
uniformly  distributed  random  variable  over  the  range  (0, 27i).  Representing  the  narrowband 

noise  w(n,  i)  by  its  complex  amplitude,  w(«) ,  we  have: 

w(n)  =  5„exp(;p„)  .  (11) 

This  is  a  complex-valued  random  variable  that  is  Gaussian-distributed  with  a  mean  that  is 


typically  zero. 

Equation  9  defines  the  complex  amplitude  of  the  signal  received  by  element  n  due 
to  the  /«th  plane  wave.  Summing  the  contributions  produced  by  a  set  of  M  plane  waves,  and 
then  adding  to  the  result  the  complex  amplitude  of  the  narrowband  noise,  we  may  express 
the  complex  amplitude  of  the  observed  signal  as: 

x(w)  =  5(«)  +  w{«) 

M 

=  Y,s(n,m)  +  w(n)  (12) 

m=\ 

M 

m=l 

We  may  put  Eq.  12  into  a  more  compact  form  by  using  matrix  notation.  First  introduce  the 
following  definitions: 


J\  n- 


N  +  l 
2  J 


+w(«),  n  =  \,2,...,N 


1.  AniVx7  vector 


16 


x= 


called  the  observed  signal  vector. 


2.  AniV^xi  vector 


^(1) 

^(2) 


x(l) 

x{2) 

x{N) 


S  = 


m=\ 

M 


m-\ 

M 


.m=l 


called  the  received  signal  vector. 


3.  An Mx 7  vector 


A  = 


ail)  ■ 

a{2) 

a(M) 


called  the  signal-in-space  vector. 


4.  An  A^xM  matrix 


(13) 


(14) 


(15) 


exp  (^jkdcosO^) 


Qxp  (^jkd cos 0  2) 


exp(7Mcos^ 


exp  [jk(N -l)d cos 0  ^)  exp  (jk(N -l)d cos 0  2)  •••  exp  (7^(iV^-l)cfcos^;j^) 

(16) 

of  propagation  delays  called  the  direction  matrix. 
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5.  AniVxl  vector 


r  w(l) 


w{N) 


called  the  noise  vector. 

Then,  using  these  definitions,  we  may  rewrite  Eq.  12  in  the  form 

X  =  S  +  W 

or,  equivalently, 

X  =  BA  +  W 


(17) 


(18) 

(19) 


where 

S  =  BA  .  (20) 

The  observed  signal  vector  X  represents  a  snapshot  of  data  corresponding  to  a  particular 
instant  of  time.  Ordinarily,  several  independent  measurements  are  made,  so  that  the  data 
available  for  processing  may  be  expressed  in  the  form: 

X(i)  =  BA(z)  +  W(/),  i  =  (21) 

where  /  is  the  total  number  of  snapshots  taken. 

Based  on  such  a  set  of  data,  we  may  use  temporal  averaging  to  improve  the 
estimation  of  the  angles  of  arrival  of  the  plane  waves  incident  upon  the  array.  The  signal-in- 
space  vector  A(0  is  modeled  as  a  stochastic  process  because  the  behavior  of  the  sources 
responsible  for  its  generation  is,  in  general,  unpredictable. 

The  spatial  correlation  matrix  of  the  observed  signal  vector  X  is  defined  by 

R  =  e[x‘(0x'(!)]  m 

where  E  is  the  expectation  operator,  the  asterisk  denotes  complex  conjugation,  and  the 
superscript  T  denotes  transposition.  Substituting  Eq.  21  in  Eq.  22,  and  recognizing  that  the 
signal-in-space  vector  A  and  the  noise  vector  W  are  statistically  independent,  we  get 
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(23) 


R  =  e[B*A*A''B'']  +  E[W*W^] 

=  B*R,B^  +  R^ 

where: 

=  spatial  correlation  matrix  of  the  signal  -  in  -  space  vector  A 

=  e[a*a^] 

R^  =  spatial  correlation  matrix  of  the  noise  vector  W 
=  E[W*W'']  . 

The  matrix  product 

R3  =  B*R,B'^  (24) 

represents  the  spatial  correlation  matrix  of  the  signal  component  of  the  signal  vector  S. 

With  the  noise  vector  W  assumed  white  noise,  we  have 

Rw  =  O-Jl  (25) 

where  cr^  is  the  variance  of  the  elemental  noise  w(n)  for  all  n,  and  I  represents  the  identity 
matrix.  Depending  on  the  structure  of  the  spatial  correlation  matrix  R^,  we  may  distinguish 

two  different  cases: 

1.  Jointly  Uncorrelated 

The  sources  responsible  for  the  signal-in-space  vector  A  are  jointly  uncorrelated  or 
noncoherent  so  that  the  off-diagonal  elements  of  the  corresponding  spatial  correlation 
matrix  are  all  zero.  Then  we  find  that  the  spatial  correlation  matrix  of  the  received  signal 
vector,  Rj  =  B*  RaB^  is  a  Toeplitz  matrix.  A  square  matrix  is  Toeplitz  if  all  the  elements  on 
its  main  diagonal  are  equal,  and  likewise  for  the  elements  on  any  other  diagonal  parallel  to 
the  main  diagonal.  The  received  signal  vector  S  is  spatially  stationary  if  its  spatial 
correlation  matrix  is  Toeplitz,  or  vice  versa.  Clearly,  if  the  spatial  correlation  matrix  R^  is 
Toeplitz,  the  spatial  correlation  matrix  of  the  observed  signal  vector,  R  is  likewise  Toeplitz. 
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2.  Correlated  Coherent 

The  sources  responsible  for  the  received  signal  vector  S  are  correlated  coherent.  The 
result  is  that  the  spatial  correlation  matrix  Rj,  is  non-Toeplitz,  and  the  received  signal  vector 
S  is  said  to  be  spatially  nonstationary.  Clearly,  If  R  is  non-Toeplitz,  the  spatial  correlation 
matrix,  R,  pertaining  to  the  observed  signal  vector,  is  also  non-Toephtz. 

C.  ARRAY  MATRIX  DECOMPOSITION  FOR  MULTIPLE  EMITTER  DOA 


For  the  location  of  the  directions  of  arrival,  the  autocorrelation  matrix  takes  the  form 
of  a  Toeplitz  matrix. 

X  =  BA  +  W 

A  =  [a(l),  a(2), . . . ,  a(  M)]^  (26) 

W  =  [  w(l),  w{2), 


Vector  X  is  explained  in  Eq.  13  and  vector  B  is  defined  in  Eq.  16;  it  contains  the 
propagation  delays,  that  is  B  =  and,  for  equally 

spaced  elements. 


~ ^Jkdicose  ■ 

1 

^jkd  cos  0 

cos  6 

^jk{N-\)d  cose 

(27) 


Thus,  the  signal  correlation  matrix  can  be  written  as 

r  =  e[bVa'b']+e[w*w'] 

=  bXb'+Rw 


(23) 


where  t  stands  for  the  complex  conjugate  transpose.  For  an  equally  spaced  array  with  N 
elements  located  at  0,  d,  2d,...,(N-l)d  the  size  of  the  correlation  matrix  in  Eq.  23  is  NxN.  In 
this  case,  the  maximum  number  of  emitters  that  can  be  detected  cannot  exceed  (N-1). 

The  matrix  R  is  Hermitian  and  always  positive  definite  if  thermal  element  noise  is 
present  Thus  R  can  be  diagonalized  using  a  rotation  matrix  E  where 
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ERE'  = 


(28) 


Xy  -  0 

•  •  • 

•  •  • 

0  •••  Xfj_ 

The  matrix  E  is  miitary  (EE'=I)  and  its  columns  are  the  eigenvectors  of  R,  where 

The  signal  covariance  matrix  is  diagonalized  by  the  same  rotation  matrix.  For  example,  if 
there  is  only  one  signal  present  (N=l),  the  weight  vector  for  the  single  signal  case  is: 

G  =  (constant)e; . 

A  measure  of  the  array  response  is  given  by  the  spectral  estimator 


m= 


1 

|g'u(0)|' 


(30) 


A  large  value  of  P  denotes  high  array  sensitivity;  that  is  a  pattern  notch  or  null.  Note 
that  this  is  not  the  same  condition  that  provides  good  signal  to  noise  for  a  conventional 
array.  If  more  than  one  emitter  is  present  the  spectral  estimator  becomes 


1 

- ^ - p  (31) 

The  power  response  for  the  conventional  equally  spaced  array  of  6  elements  using  ordinary 
beamforming  is  shown  in  Figure  7.  In  this  case  there  is  no  ability  to  resolve  two  signals,  the 
null  locations  are  fixed  relative  to  each  other.  Figure  8  shows  the  response  of  a 
Caratheodory  array  [Ref.  1]  with  four  elements  located  at  0,  d,  5dl2,  and  'id.  The  emitter 
directions  are  at  20,  50,  90,  130  and  160  degrees  respectively.  Note  that  all  of  the  emitters 
are  essentially  resolved,  but  the  sensitivity  is  only  lOdB  in  some  cases. 
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Figure  7.  Conventional  Equally  Spaced  Array  with  Six  Elements  when  the  Emitter 

Directions  are  20, 50, 90, 130,  and  160 
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Figure  8.  Conventional  Processing  Scheme  for  a  Caratheodory-array  with  Four 
Elements  when  the  Emitter  Directions  are  20, 50, 90, 130,  and  160 
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IV.  RESIDUE  ANTENNA  ARCHITECTURE 


A.  OVERVIEW  OF  RESIDUE  NUMBER  SYSTEMS 

In  the  first  century  A.D.  the  Chinese  scholar  Sun  Tzu  authored  a  book  posing  the 
problem  of  determining  a  number  having  remainder  2, 3,  and  3  when  divided  by  3, 4,  and  5 
respectively.  The  answer  is  23  and  is  found  by  using  a  three-moduli  (3,4,5)  number  system. 
Later  this  theorem  was  called  the  Chinese  Remainder  Theorem  [Ref  8].  Eventually  this 
theorem  became  an  important  cornerstone  in  the  modem  theory  of  residue  number  system 
(RNS)  arithmetic. 

Digital  systems  that  use  residue  arithmetic  units  play  an  important  role  in  ultra-high 
speed,  dedicated,  real-time  systems  that  support  pure  parallel  processing  of  integer  valued 
data  Residue  arithmetic  performs  addition,  subtraction,  and  multiplication  as  concurrent 
operations,  side-stepping  one  of  the  principal  arithmetic  delays:  managing  carry 
information.  Some  applications  such  as  the  fast  Fourier  transform,  (iigital  filtering,  and 
image  processing  utilize  the  efficiencies  of  RNS  arithmetic  in  addition  and  multiplication. 

The  RNS  has  been  employed  efficiently  in  the  implementation  of  digital  processors. 
The  RNS  is  composed  of  a  number  of  moduli  m-^.  The  integers  within  each  RNS  modulus 
are  representative  of  a  sawtooth  waveform  with  the  period  of  the  waveform  equal  to  the 
modulus. 

Figure  9  shows  the  RNS  folding  waveform  and  the  corresponding  RNS  states  for 
the  moduli  /m,=3,  nif^,  mf=5.  In  this  scheme,  given  the  integer  values  within  the 
individual  modulus  are  given  by  the  row  vector: 

-l]  (32) 

The  required  number  of  comparators  for  the  zth  channel  is  m-l .  Notice  that  all  the  integers 
within  Eq.  32  form  a  complete  system  of  length  Wj .  The  dynamic  range  of  this  scheme,  as 
well  as  the  number  of  quantization  levels,  is  given  by: 
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(33) 


N 

Dynamic  Range  =  n  m,. 

j=i 


Table  1  illustrates  the  RNS  encoding  values  for  the  moduli  (3,4,5). 


^1=3 


«*2=4 


wz3=5 


Figure  9.  RNS  Folding  Waveforms  and  Output  Code  for  the  Moduli  (3,4,5) 
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Analog  Moduli 

Input  »i,=3  mi=4  fM3=5 
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1 

4 

50 

2 

2 

0 

51 

0 

3 

1 

52 

1 

0 

2 

53 

2 

1 

3 

54 

0 

2 

4 

55 

1 

3 

0 

56 

2 

0 

1 

57 

0 

1 

2 

58 

1 

2 

3 

59 

2 

3 

4 

60 

0 

0 

0 

Table  1.  RNS  Encoding  Procedure  for  the  Moduli  (3,4>5) 


B.  RNS  ANTENNA  ARCHITECTURE 

In  the  RNS  antenna  architecture  shown  in  Figure  10,  phase  detectors  are  used  to 
combine  the  individual  pairs  of  elements.  The  comparator  threshold  levels  and  subsequent 
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logic  block  are  particular  to  the  RNS.  The  phase  response  (A^)  of  each  interferometer 
corresponds  to  a  sawtooth  folding  waveform. 


Figure  10.  RNS  Antenna  Architecture 
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For  a  given  modulus  nip  the  number  of  folds  needed  to  apply  the  RNS  code  within 
the  dynamic  range  is  given  by: 

dynamic  range  M 


”i,RNS  ~ 


m.. 


m. 


(34) 


From  Figure  10 


A^  =  0°  , 

(9  =  90°=>  = — 


(35) 


0  = -90°=>  A^zi  = - 


Ind 

~Y 


For  a  DOA  of  180°,  the  phase  response  corresponds  to  a  change  in  phase  of  271.  We  require 
that 

A<^90°)  -  A<i5(-90'’)  =  Ik  ,  (36) 

which  implies 


Aiid  ^  ,  'k 

- =  2%,  or  d  =  —  . 

X  2 

Finally,  the  distance  between  elements  is  given  by: 

^X^  MfV 

“/,RNS  “  ”i,RNS 


(37) 


(38) 


\2J  m,\2J' 

Using  the  same  moduli  set  with  center  frequency  of  /o  =  8.5  GHz ,  the  distances  for  each 


pair  of  elements  are: 


Ji=10X,  , 

=  7.5A,  , 

=  6.0A,  . 


These  interferometer  distances  are  used  to  investigate  numerically  the  performance  of  the 
RNS  antenna  architecture. 
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C.  RNS  ARRAY  CORRELATION  MATRIX  FOR  MULTIPLE  EMITTERS 

In  this  section  we  briefly  investigate  the  capability  of  the  RNS  array  with  /n,=3, 
W2=4,  and  mf^5  to  resolve  multiple  emitters  at  the  same  frequency.  The  same  approach  is 
used  for  as  the  minimum  redundancy  array  (Chapter  2,  Section  C).  For  a  mmimum 
redundancy  array,  not  all  elements  must  be  present,  but  all  correlation  lags  are  present  (d, 
2d,  3d„  ...,(N-l)d)  since  all  spacings  are  present  at  least  once.  Therefore  the  array  correlation 

matrix  can  be  filled  and  consequently  up  to  — ^  emitters  can  be  resolved  with  N 
elements. 

For  the  RNS  array  not  all  elements  spacings  are  present.  In  the  case  of  mi=3,  W2=4 
and  mf=5  the  spzicings  that  are  present  are  Od,  \2d  (6 A),  \5d  (7.5 A),  20c?  (10 A),  3d 
(1 .5  A ),  5c?  (2.5  A  ),  and  8c?  (4  A )  as  shown  in  Figure  1 1 . 


10A,(20d) 


- 7.51(15d) 

62(12d) - 


V  V  V 


r  1 

X  >1 

< - 

^2.5 

-4  A. - > 

Figure  11.  RNS  Antenna  Spacing  for  the  (3,4,5)  Array 
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The  standard  periodic  array  would  have  a  completely  filled  21  x  21  correlation 
matrix.  Not  all  spacings  are  represented  in  the  RNS  array,  and  if  the  missing  correlation  lags 
are  set  to  zero  (or  noise)  the  matrix  would  be: 


R  = 


Ro  0  0  R3  0  Rs  0  0  Rg  0 

0  Ro  0  0  R3  0  R5  0  0  Rg 

0  0  Ro  0  0  Rj  0  R5  0  0 

R*3  0  0  Ro  0  0  R3  0  Rs  0 

0  R*3  0  0  Ro  0  0  R3  0  Rj 

R*5  0  R*3  0  0  Ro  0  0  R3  0 

0  R*5  0  R’3  0  0  Ro  0  0  R3 

0  0  R*5  0  R*3  0  0  Ro  0  0 

R*8  0  0  R*5  0  R*3  0  0  Ro  0 

0  R*8  0  0  R*5  0  R*3  0  0  Ro 

0  0  R*8  0  0  R*5  0  R*3  0  0 

0  0  0  R*8  0  0  R*5  0  R*3  0 

R*i2  0  0  0  R'g  0  0  R*5  0  R*3 

0  R*i2  0  0  0  R*8  0  0  R*5  0 

0  0  R*,2  0  0  0  R'g  0  0  R'5 

R',5  0  0  R*i2  0  0  0  R'g  0  0 

0  R*  5  0  0  R',2  0  0  0  R'g  0 

0  0  R'i5  0  0  R'i2  0  0  0  R'g 

0  0  0  R',5  0  0  R*,2  0  0  0 

0  0  0  0  R'i5  0  0  R'i2  0  0 

R*2o  0  0  0  0  R'k  0  0  R'i2  0 


0  0  Ri2  0  0  R,5  0  0  0 

0  0  0  R,2  0  0  R,5  0  0 

Rg  0  0  0  R,2  0  0  R,5  0 

0  Rg  0  0  0  Ri2  ^  ^  Ri5 

0  0  Rg  0  0  0  Ri2  0  0 

R5  0  0  Rg  0  0  0  Ri2  0 

0  R5  0  0  Rg  0  0  0  Ri2 

R3  0  Rs  0  0  Rg  0  0  0 

0  R3  0  R5  0  0  Rg  0  0 

0  0  R3  0  R5  0  0  Rg  0 

Ro  0  0  R3  0  R5  0  0  Rg 

ORoO  OR3OR5O  0 
0  ORoO  OR3OR5O 
R'3  0  0  Ro  0  0  R3  0  R5 

0  R'3  0  0  Ro  0  0  R3  0 

R's  0  R'3  0  0  Ro  0  0  R3 

0  R'5  0  R'3  0  0  Ro  0  0 

0  0  R'5  0  R'3  0  0  Ro  0 

R'g  0  0  R*5  0  R*3  0  0  Ro 

0  R'g  0  0  R'5  0  R'3  0  0 

0  0  R'g  0  0  R'5  0  R‘3  0 


0  R20 
0  0 
0  0 
0  0 
R,5  0 

0  R,5 
0  0 
Ri2  0 
0  Ri2 
0  0 
0  0 
Rs  0 
0  Rg 
0  0 
R5  0 
0  R5 
R3  0 
0  R3 
0  0 
Ro  0 
0  Ro 


Figure  12  shows  the  array  response  for  50,  90  and  130  degrees.  As  expected  the 
performance  is  not  good  because  of  the  large  number  of  zeros  in  the  array . 
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Figure  12.  RNS  Array  Response  for  Multiple  Emitters  when  the  Emitter  Directions 

are  50, 90,  and  130 

D.  COMPUTER  SIMULATION  OF  RNS  CODE 

The  quantization  of  the  phase  response  from  each  pair  of  elements  is  made  separately  using 
a  small  bank  of  comparators.  The  threshold  levels  for  each  comparator  are  shown  in  Table 
2.  If  the  phase  response  is  plotted  as  a  function  of  sin0,  an  ideal  sawtooth  waveform  is 
obtained  as  shown  at  the  top  of  Figure  13.  An  actual  plot  of  the  phase  response  is  shown  in 
Figure  14  as  a  function  of  the  input  direction  of  arrival  0.  The  sawtooth  waveform  stretches 
out  as  0,  the  input  direction  of  arrival,  goes  from  0  to  90°.  For  a  specific  input  direction 
of  arrival,  there  will  be  a  corresponding  integer  number  of  comparators  in  the  “on” 
state  within  the  modulus  w,.  The  value  of  the  integers  depends  on  the  measured  phase 
difference  associated  with  each  pair  of  elements  for  the  specific  input  direction  of  arrival. 
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The  three  RNS  codes  (one  for  each  pair  of  antenna  elements)  are  indices  to  an  RNS  logic 
look-up  table  that  gives  the  resolved  direction  of  arrival.  Table  3  shows  the  three  RNS 
codes  and  the  resolved  direction  of  arrival.  The  numbered  codes  represent  the  number  of 
comparators  in  the  “on”  state. 


Moduli 

T, 

Tz 

T, 

T4 

3 

-1.047 

1.047 

4 

-1.571 

0.00 

1.571 

5 

-1.885 

-0.628 

0.628 

1.885 

Table  2.  Threshold  Levels  for  the  Moduli  (3,4,5) 


Input  Direction  of  Arrival  (Sin  0) 

Figure  13.  Phase  Response  as  a  Function  of  sin0  for  the  Moduli  (3,4,5) 
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Figure  14.  Phase  Response  as  a  Function  of  the  Input  Direction  of  Arrival  (0) 

for  the  Moduli  (3,4)5) 

Figure  15,  shows  a  plot  of  the  resolved  direction  of  arrival  as  a  function  of  0,  the 
input  direction  of  arrival.  The  RNS  antenna  system  gives  the  direction  of  arrival 
information  over  the  complete  field  of  view  without  ambiguity  (DOA  increment  is  0.3°). 
However  the  accuracy  varies  depending  on  the  actual  direction  of  arrival  of  the  incoming 
signal.  Notice  also  that  the  quantization  width  increases  as  the  direction  of  arrival 
approaches  ±90°. 


32 


33 


Resolved  Direction  of  Airivai  ( degrees ) 


Figure  15.  Resolved  Direction  of  Arrival  in  Function  of  the  Input  Dire< 

when  DOA  Increment  is  0.3° 
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V.  ENCODING  ERROR  CORRECTION 


A.  ENCODING  ERRORS 

The  DOA  estimates  previously  shown  in  Figure  15  are  based  on  ideal  hardware 
performance.  One  significant  assumption  is  that  all  of  the  comparator  thresholds  are  ideal 
and  cross  the  felling  phase  response  waveforms  together  at  each  DOA  code  transition  point 
If  the  samples  are  generated  by  stepping  the  emitter  direction  between  -90°  and  90°  in 
increments  of  0.03°,  for  example,  some  samples  will  fall  about  the  code  transition  point  and 
result  in  encoding  errors  due  to  some  comparators  do  not  switch  simultaneously.  When  the 
emitter  is  swept  through  an  angular  sector,  the  DOA  estimates  from  a  realizable  system 
could  contain  large  spikes  or  “glitches”  at  these  points.  An  illustration  of  the  encoding 
errors  is  given  in  Table  4.  The  columns  in  Table  4  correspond  to  the  first  5  RNS  codes 
shown  in  Table  3  with  each  code  sampled  three  times.  The  seventh  sample  shows  the  nature 
of  an  encoding  error  where  some  of  the  comparators  in  modulus  5  do  not  switch  at  the  same 
time.  The  seventh  sample  should  be  mfO.,  mi=2  and  mi=2.  This  results  in  a  glitch  at  this 
sample  within  the  transfer  function  as  shown  in  Figure  16. 


Sample 

1 

m 

m 

5 

6 

:  7 

o 

10 

11 

m 

13 

14 

Modulus  3 

0 

o 

o 

a 

1 

1 

2 

B 

B 

B 

0 

B 
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1 

Modulus  4 

O 

0 

o 

a 

1 

1 

2 

B 

2 

B 

3 

B 

0 

0 

Modulus  5 

O 

o 

_L 

D 

1 

1 

B 

2 

B 

3 

B 

:  4 

4 

Table  4.  Encoding  Error  in  the  RNS  Output  Code 


To  examine  the  nature  of  the  encoding  errors,  a  simulation  was  performed  with  the 
threshold  levels  shifted  from  the  ideal  (error-free)  threshold  levels  in  Table  2  to  a  set  of  new 
threshold  levels.  The  errors  are  induced  in  order  to  compare  three  interpolation  methods 
designed  to  eliminate  the  encoding  errors.  The  new  threshold  values  are  shown  in  Table  5, 
and  the  glitches  in  the  DOA  are  shown  in  Figure  17,  with  the  total  number  of  glitches  equal 
to  44. 
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Figure  16.  RNS-to-Decimal  Logic  Block  and  Glitch 
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Moduli 

T, 

Tj 

T3 

T4 

3 

-1.047 

1.047 

4 

-1.560 

0.000 

1.560 

5 

-1.870 

-0.640 

0.640 

1.870 

Table  5.  New  threshold  Levels  for  the  Moduli  (3,4,5) 


Encoding  Error  in  Quatization  (glitch=  44, sampling  distance:ds=0.03  degree) 
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Figure  17.  Glitches  and  Resolved  Direction  of  Arrival 


B.  ADDITIONAL  COMPARATORS  TO  ISOLATE  ENCODING  ERRORS 

The  original  threshold  levels  within  the  minimum  modulus  are  set  to  the  level  of 
1.047,  -1.047  radians.  The  phase  response  for  the  minimum  modulus  (m,=3)  is  shown  in 
Figure  1 8  and  illustrates  the  position  of  the  comparator  threshold  levels. 
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Figure  18.  Original  Threshold  Level  in  Modulus  3 

To  isolate  the  possible  encoding  errors,  2w,  comparators  are  used  in  the  minimum 
modulus  (instead  of  -1).  The  comparator  thresholds  levels  are  such  that  a  small  band  is 
set  up  about  the  code  transition  points.  The  values  of  the  threshold  levels  depend  on  the 
width  of  the  band  and  can  be  expressed  ais  a  percentage  of  the  least  significant  angle  (LSA) 
width.  The  six-threshold  levels  for  modulus  mf^'i  are  shown  in  Figure  19,  and  listed  in 
Table  6  as  a  function  of  the  LSA  percentage. 

The  encoding  errors  are  isolated  by  comparing  the  parity  value  for  the  2m^ 
comparators  in  the  minimum  modulus.  If  the  parity  value  is  odd  number  (1,  3  and  5),  then 
the  sample  falls  outside  the  small  transition  band  and  does  not  present  a  problem.  If  the 
parity  is  an  even  number  (2, 4  and  6),  the  sample  represents  a  possible  DOA  encoding  error. 
The  goal  of  the  interpolation  method  is  to  determine  the  best  DOA  estimate  given  the  fact 
that  die  sample  could  present  a  possible  error. 
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Plase  Output 


Figure  19.  Additional  Threshold  Levels  and  RNS  Output  Code 


LSA  Percentage 

T„ 

Tas 

T.4 

T.S 

T.6 

10.00 

-3.0369 

-1.1519 

-0,9425 

0.9425 

1.1519 

3.0369 

9.00 

-3.0473 

-1.1414 

-0.9529 

0.9529 

1.1414 

3.0473 

8.00 

-3.0578 

-1.1310 

-0.9634 

0.9634 

1.1310 

3.0578 

7.00 

-3.0683 

-1.1205 

-0.9739 

0.9739 

1.1205 

3.0683 

6.00 

-3.0788 

-1.1100 

-0.9844 

0.9844 

1.1100 

3.0788 

5.00 

-3.0892 

-1.0996 

-0.9948 

0.9948 

1.0996 

3.0892 

4.00 

-3.0997 

-1.0891 

-1.0053 

1.0053 

1.0891 

3.0997 

3.00 

-3.1102 

-1.0786 

1.0786 

3.1102 

2.00 

-3.1206 

-1.0681 

-1.0263 

1.0263 

1.0681 

3.1206 

1.00 

-3.1311 

-1.0577 

-1.0367 

1.0367 

1.0577 

3.1311 

Table  6.  LSA  Percentage  and  Six- Additional  Threshold  Values 
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C  EVTERPOLATION 


A  block  Hiflgram  of  the  RNS  antenna  architecture  including  the  interpolation 
processing  is  shown  in  Figure  20. 
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Direction  of  Arrival  9 


Figure  20.The  Procedure  of  DOA  Decoder 

The  RNS-to-binary  logic  block  is  set  after  the  parity  circuit  as  the  interpolation 
procedures  are  accomplished  in  the  binary  number  system  as  described  below.  The  RNS 
code  and  corresponding  binary  values  are  shown  in  Table  7.  The  encoding  errors  at  samples 
7  and  13,  and  the  RNS  code  in  Table  4  are  converted  to  a  six-bit  binary  code  and  are  shown 
in  Table  8. 
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Angle 

RNS 

Code 

Six-Bit 

Binary 

Code 

Angle 

RNS 

Code 

Six-Bit 

Binary 

Code 

-82.5824 

0 

0 

0 

000000 

0.9551 

0 

2 

0 

011110 

-72.0627 

1 

1 

1 

000001 

2.8664 

1 

3 

1 

011111 

-66.5593 

2 

2 

2 

000010 

4.7809 

2 

0 

2 

100000 

-62.1158 

0 

3 

3 

000011 

6.7007 

0 

1 

3 

100001 

-58.2581 

1 

0 

4 

000100 

8.6382 

1 

2 

4 

100010 

-54.7864 

2 

1 

0 

000101 

10.5655 

2 

3 

0 

100011 

-51.5928 

0 

2 

1 

000110 

12.5152 

0 

0 

1 

100100 

-48.6111 

1 

3 

2 

000111 

14.4797 

i 

1 

2 

100101 

-45.7968 

2 

0 

3 

001000 

16.4618 

2 

2 

3 

100110 

-43.1187 

0 

1 

4 

001001 

18.4644 

0 

3 

4 

100111 

-40.5534 

1 

2 

0 

001010 

20.4907 

i 

0 

0 

101000 

-38.0832 

2 

3 

1 

001011 

22.5442 

2 

1 

1 

101001 

-35.6940 

0 

0 

2 

001100 

24.6287 

0 

2 

2 

101010 

-33.3745 

1 

1 

3 

001101 

26.7487 

i 

3 

3 

101011 

-31.1155 

2 

2 

4 

001110 

28.9091 

2 

0 

4 

101100 

-28.9091 

0 

3 

0 

001111 

31.1155 

0 

1 

0 

101101 

-26.7487 

1 

0 

1 

010000 

33.3745 

1 

2 

1 

101110 

-24.6287 

2 

1 

2 

010001 

35.6940 

2 

3 

2 

101111 

-22.5442 

0 

2 

3 

010010 

38.0832 

0 

0 

3 

110000 

-20.4907 

1 

3 

4 

010011 

40.5534 

1 

1 

4 

110001 

-18.4644 

2 

0 

0 

010100 

43.1187 

2 

2 

0 

110010 

-16.4618 

0 

1 

1 

010101 

45.7968 

0 

3 

1 

110011 

-14.4797 

1 

2 

2 

010110 

48.6111 

1 

0 

2 

110100 

-12.5152 

2 

3 

3 

010111 

51.5928 

2 

1 

3 

110101 

-10.5655 

0 

0 

4 

011000 

54.7864 

0 

2 

4 

110110 

-8.6382 

1 

1 

0 

oiioof 

58.2581 

1 

3 

0 

110111 

-6.7007 

2 

2 

1 

IB 

0 

1 

111000 

-4.7809 

0 

3 

2 

011011 

66.5593 

0 

1 

2 

111001 

-2.8664 

1 

0 

3 

011100 

72.0627 

1 

2 

3 

111010 

-0.9551 

2 

1 

4 

011101 

82.5824 

iKKi 

111011 

Table  7.  RNS  code  and  Six-bit  Binary  Code  in  Moduli  (3,4,5) 
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Sample 

1 

2 

3 

4 

5 

6 

7 

9 

10 

LiL 

la 

14 

MSB 
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0 

0 

0 

0 

0 

0 

0 
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■■ 
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■SB': 

0 

- 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

- 

0 

0 

0 

0 

0 

0 

0 

11 

0 

0 

0 

0 

1 

0 

- 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

,  1 

1 

1 

1 

n 

0 

0 

LSB 

0 

0 

0 

cu 

m 

□2 

mm 

11 

0 

1 

1 

0 

0 

MSB  represents  the  Most  Significant  Bit  of  the  converter 


LSB  is  the  Least  Significant  Bit  in  the  six-bit  converter 

Table  8.  Encoding  Error  in  Six-bit  Binary  Code 

1.  LSB  Shift  Method 

In  Table  4  glitches  occured  at  samples  7  and  13.  From  that  binary  data  the  output 
signal  jumps  from  the  value  given  by  000001  through  000110  to  000010  as  shown  in  Table 
8.  To  resolve  these  glitches,  we  investigate  the  use  of  a  7x5  Error  Correction  Block  (ECB) 
matrix.  The  output  of  a  parity  circuit  determines  whether  the  sample  is  good  or  bad.  If  the 
parity  output  is  an  odd  number  (1,  3  and  5)  the  sample  is  good,  whereas  an  even  number  (2, 
4  and  6)  the  output  represents  a  possible  bad  sample.  The  parity  circuit  and  error  correction 
block  are  shown  in  Figure  21 . 

a.  LSB  Shift  Algorithm 

The  parity  and  the  6-bit  word  (from  the  PLA)  are  clocked  in  parallel  through 
a  5-state  buffer.  The  buffer  allows  the  6-bit  word  to  be  replaced  with  an  interpolated  value 
depending  on  the  value  of  the  parity.  If  the  parity  is  odd  (good  sample),  the  6-bit  word  is 
clocked  through  the  buffer  unaltered.  If  the  parity  is  even,  the  values  ahead  and  behind  the 
even  parity  value  are  used  to  interpolate  the  possible  encoding  error.  The  interpolation 
process  compares  the  most  significant  bits  (MSBs)  of  the  odd  parity  words  stopping  at  the 
point  where  the  MSBs  do  not  agree.  The  MSBs  that  are  the  same  are  copied  into  the  MSBs 
for  the  even  parity  word.  The  LSBs  come  directly  from  the  original  (even  parity)  word. 
Figure  22  shows  an  example  where  there  is  one  out  of  five  possible  bad  samples.  In  the 
interpolation  process,  P2  and  Pq  combine  to  form  the  MSBs  of  T]  (an  intermediate  state). 
The  LSBs  of  T,  come  from  Pj.  At  time  t+1,  T,  is  shifted  into  Pj.  The  situation  where  two 
samples  are  even  parity  is  shown  in  Figure  23.  In  this  case  P3  and  Pq  combine  to  form  the 
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MSBs  of  Tj  and  T2.  The  LSBs  of  T,  and  Tj  are  copied  directly  from  P,  and  Pj.  Tj  is  shifted 
into  P2  and  T2  is  shifted  into  P3  at  t+1.  The  interpolation  procedure  for  three  out  of  five 
samples  having  even  parity  is  shown  in  Figure  24.  The  interpolation  procedure  for  four  out 
of  five  samples  having  even  parity  is  shown  in  Figure  25. 


Parity  =  1  (Sample  OK) 
Output  0  (Sample  Bad) 


7x5  Error  Correction  Block 


> 


LSB _ 

Sn=  Shift  Register 


Figure  21.  Parity  Circuit  and  Error  Correction  Block 
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Figure  22.  LSB-Shift  Method  in  Case  of  1/5  Error  Samples 


in  Case  of  3/5  Error  Samples 


b.  Simulation  Results 

The  DOA  transfer  function  including  the  encoding  error;  was  shown  in 
Figure  17.  The  LSB  shift  method  with  the  interpolation  width  equal  to  1%  of  the  LSA  is 
shown  in  Figure  26.  The  total  number  of  glitches  has  been  reduced  from  44  to  42.  Figure  27 
through  Figure  30  show  the  result  of  the  LSB  shift  method  with  the  interpolation  width 
ranging  from  2%  to  5%  respectively.  Figure  31  summarizes  the  results  and  shows  the  total 
number  of  glitches  as  a  function  of  the  isolation  bandwidth  (percentage  of  LSA).  The  larger 
the  bandwidth,  the  fewer  the  number  of  glitches. 


Interpolation  Resutti  %  LSB’Shift  Method  ( glitch=  42,  sampling  distancerds^O.OS  degree) 


-80  -60  -40  -20  0  20  40  60  80 

Input  Direction  of  Arrival  ( degrees ) 


Figure  26.  LSB-Shift  Interpolation  Results,  1%  of  LSA 
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Figure  28.  LSB-Shift  Interpolation  Results,  3%  of  LSA 
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Revived  pOA  and  Quantization  Error  ( degrees )  Resolved  DOA  and  Quantization  Error  ( degrees ) 


Interpolation  Result, 4  %  LSB-ShifI  Method  ( glitch=  15,  sampling  di8tance:ds=0.03  degree 


Figure  29.  LSB-Shift  Interpolation  Results,  4%  of  LSA 


Interpolation  Results  %  LSB-Shift  Method  ( glitch=  8,  sampling  dlstance:ds=0.03  degree) 


Figure  30.  LSB-Shilt  Interpolation  Results,  5%  of  LSA 
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Figure  31.  Number  of  Glitches  vs  LSA  Percentage  After  LSB-Shift  Method 

2.  Random  LSB  Method 

The  random  LSB  method  compares  the  most  significant  bits  (MSBs)  of  the  odd 
parity  words  stopping  at  the  point  where  the  MSBs  do  not  agree.  The  MSBs  that  are  the 
same  are  copied  into  the  MSBs  for  the  even  parity  word.  The  LSBs  do  not  come  fi'om  the 
original  (even  parity)  word.  Instead  the  LSBs  are  chosen  randomly  (0  or  1).  The 
interpolation  procedure  for  three  out  of  five  samples  having  even  parity  is  shown  in  Figure 
32. 

The  simulation  for  the  random  LSB  method  with  the  interpolation  width  changing 
from  1%  to  5%  are  shown  in  Figure  33  to  Figure  37.  The  total  number  of  glitches  changed 
randomly,  so  it  is  difficult  to  determine  any  relationship  between  the  LSA  percentage  and 
glitches.  The  results  are  summarized  in  Figure  38  and  the  randomizing  of  the  LSBs  seems 
to  create  more  problems. 
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Figure  32.  Procedure  of  Random  LSB  Method 


Interpolation  Result, 1  %  Random-LSB  Method  ( glitch=  82,  sampling  distance:ds=0.03  degrees 


Interpolation  Result, 2  %  Random-LSB  Method  { glltch=  95,  sampling  distance:ds=0.03  degrees 


Figure  34.  Random-LSB  Method  Result,  2%  of  LSA 
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Resolved  DOA  and  Quantization  Error  ( degrees ) 


Interpolation  Result, 3  %  Random-LSB  Method  ( glitch=  1 14.  sampling  distancetdssO.OS  degree) 


Figure  35.  Random-LSB  Method  Result,  3%  of  LSA 

InterpolaUon  Result, 4  %  Random-LSB  Method(  glitch*  sa  sampling  distanco:ds=0.03  degree) 


Figure  36.  Random-LSB  Method  Result,  4%  of  LSA 


Number  of  Cliches 


Irrterpoiation  Result, 5  %  Random4.SB  Method  ( g]itch=  67.  sampling  dt8tanc8:ds=0.03  degree) 


Figure  37.  Random-LSB  Method  Result,  5%  of  LSA 


Figure  38.  Number  of  Glitches  vs  LSA  Percentage  after  Random  LSB  Method 


3.  Shift  Last  Good-Sample  Method 

In  the  last  two  methods,  some  error  samples  cannot  be  resolved  completely.  This  is 
because  some  portion  of  least  significant  bits  remains  in  the  error  correction  block  buffer 
without  changing.  With  the  previous  methods  some  glitches  may  remain  in  the  sample. 
However,  the  shift  last  good-sample  method  removes  the  encoding  error  perfectly  after 
interpolation  by  shifting  directly  fi-om  good  sample's  6-bit  value  to  error  sample's  6-bit 
value.  The  procedure  of  this  method  is  shown  in  Figure  39. 
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Figure  39.  Procedure  of  Shift  Last  Good-Sample  Method 


Interpolation  Result, 1  %  Shift  Last-Good  Bit  Method  (glitch=  36,  sampling  distance;ds=0.03  degrees ) 


Figure  40.  Shift  Last  Good-Sample  Method  Result,  1%  of  LSA 


Interpolation  Result, 2  %  Shift  Last-Good  Bit  Method  (glitch=  23,  sampling  dlstance:ds=0.03  degrees ) 


Figure  41.  Shift  Last  Good-Sample  Method  Result,  2%  of  LSA 
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Interpolation  Re$ult,3  %  Shift  Last<Cood  BH  Method  (glitch^  9,  sampling  distance:d$»0.03  degrees ) 


Figure  42.  Shift  Last  Good-Sample  Method  Result,  3%  of  LSA 

Interpolation  R0sult,4  %  Shift  Last-Good  Bit  Method  (glitch=  7,  sampling  distanc8:ds=0.03  degrees ) 


Figure  43.  Shift  Last  Good-Sample  Method  Result,  4%  of  LSA 
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Interpolation  Results  %  Shift  Last-Good  Bit  Method  (glitch=  3,  sampling  di8tancG:ds=0.03de^ees ) 


Figure  44.  Shift  Last  Good-Sample  Method  Result,  5%  of  LSA 


Percentage  of  LSA 


Figure  45.  Number  of  Glitches  vs  LSA  Percentage  After  Shift  Last  Good-Sample 

Method 
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VI.  CONCLUDING  REMARKS 


The  important  contribution  of  this  thesis  is  detailing  a  residue  antenna  architecture 
and  interpolation  methods  that  may  hold  the  key  to  resolving  the  encoding  errors  (glitches) 
using  an  error  correction  block.  The  essential  part  of  interpolation  procedure  is  to  determine 
the  additional  threshold  levels  according  to  LSA  percentage  in  the  first  modulus  (mi=3),  so 
that  the  parity  circuit  can  be  applied  to  the  error  correction  logic  block  system.  Using  the 
residue  number  system  and  the  six-bit  binary  code  system  it  is  very  fast  and  develops  good 
resolution  for  converting  the  analog  input  into  the  direction  of  arrival  (DO  A)  output. 

The  three  mterpolation  method  results  show  that  the  number  of  glitches  is 
determined  by  the  percentage  of  LSA  in  three  trials.  In  the  LSB  shift  method  and  shift  last 
good  sample  method,  the  simulation  results  indicate  that  the  higher  the  LSA  percentage  the 
better  the  result  (i.e.  fewer  glitches  after  interpolation).  The  random  LSB  method  shows  no 
trends  between  LSA  percentage  and  the  interpolation  results  because  of  the  number  of 
glitches  is  randomly  changed.  A  comparison  of  three  interpolation  methods  is  shown  in 
Figure  46. 

The  simulated  results  indicate  that  the  RNS  provides  a  more  robust  solution  to  high- 
resolution  DF  and  also  has  excellent  performance.  In  fact,  since  the  phase  response  has  a 
"uniform"  distribution  as  a  function  of  angle  of  arrival  it  is  possible  to  obtain  the  direction 
of  arrival  information  over  the  entire  field  of  view.  To  avoid  ambiguities  in  the  resolved 
direction  of  arrival  the  spacing  between  elements  must  be  established  for  the  highest 
firequency  emitter.  On  the  other  hand,  the  frequency  of  the  incoming  signal  has  to  be  known 
in  order  to  make  the  necessary  corrections  in  the  RNS  logic  block  using  a  fast  correction 
algorithm  to  obtain  the  correct  direction  of  arrival.  The  performance  of  the  algorithm  was 
quantified  and  shown  to  have  excellent  results  over  a  large  bandwidth.  Future  work  should 
include  a  demonstration  array  with  both  the  microwave  and  A/D  hardware  incorporated  as 
shown  in  Figure  10. 
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Percentage  of  LSA 

Figure  46.  Comparison  of  the  Three  Interpolation  Method  Results 


APPENDIX  MATLAB  CODE 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  This  program  'Phase3'  calculates  the  Direction  Of  Arrival  (DO A)  using  % 

%  the  phase  response  of  the  array  and  RNS  code  in  fast  A/D  converter  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  program  1 
clear  all 

%  Calculate  the  Dynamic  Range 

ml=3; 

m2=4; 

m3=5; 

M=ml*m2*m3; 

rad=pi/180; 

ds=3\input('What  is  the  step  size  (degree)  ? ') 

LSA_percent  =input('What  is  the  percentage  of  LSA  ?(%)') 
fo=input('What  was  the  frequency  of  the  antenna  (GHz)  ?') 
fr==input(''What  is  the  frequency  of  the  received  signal(GHz)?') 
fo=fo*le9; 
fr=fr*le9; 

theta=linspace(-90,90,M/ds); 
thet=sin(theta*rad)*  90; 
c=3e8; 

L=c/fo; 

Ll=c/fr; 

%  Distances  Between  Each  Pair  of  Elements 

dl=0.5*M/ml*L; 

d2=0.5*M/m2*L; 

d3=0.5*M/m3*L; 

%  Generates  the  Phase  Response  of  Each  Pair  of  Elements 
phi  1  =(2*pi*dl/Ll  )*sin(theta*rad); 
modlp=angle(exp(j*((phil+pi)))); 

phi2=(2*pi*d2/Ll)*sin(theta*rad); 
mod2p=angle(exp(j  *((phi2)))); 
phi3=(2*pi*d3/Ll)*sin(theta*rad); 
mod3p=angle(exp(i*((phi3)+pi))); 

%  Find  the  Threshold  Values  for  Each  Modulus 
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thresl 


%  Generates  the  Comparators  OuQ)ut  RNS  code  for  Modulus  ml=3,m2=4,m3=4 

for  i=l  ;M/ds 

ifmodlp(i)<=Tl(l) 

modl(i)=0; 

elseif  modlp(i)>Tl(2) 

modl(i)=2; 

else 

modl(i)=l; 

end 

ifmod2p(i)<=T2(l) 

mod2(i^0; 

elseif  mod2p(i)<=T2(2) 
mod2(i)=l; 

elseif  mod2p(i)<=T2(3) 

mod2(i)=2; 

else 

mod2(i)=3; 

end 

ifmod3p(i)<=T3(l) 

mod3(i)=0; 

elseif  mod3p(i)<=T3(2) 
mod3(i)=l; 

elseif  mod3p(i)<=T3(3) 
mod3(i)=2; 

elseif  mod3p(i)<=T3(4) 

mod3(i)=3; 

else 

mod3(i)=4; 

end 

end 

figure(l) 
whitebg 
subplot(61 1) 
plot(theta,mod  1  p) 
axis([-90,90,-4,4]); 

title(['Foldingwave  and  RNS  output  code  (sampling 

distance:ds-,num2str(ds*3), ’degree)']),',  num2str(l/ds),'samples  in  3  degree )']) 
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grid 

subplot(612) 
plot(theta,mod2p) 
axis([-90,90,-4,4]); 
ylabel('Folding  waves ') 
grid 

subplot(613) 

plot(theta,mod3p) 

axis([-90,90,-4,4]); 

grid 

subplot(614) 
stairs(theta,modl ) 
axis([-90,90,-0.2,2.2]); 
grid 

subplot(615) 
stairs(theta,mod2) 
axis([-90,90,-0.2,3-2]); 
ylabel('RNS  Output  Code') 
grid 

subplot(616) 

stairs(theta,mod3) 

axis([-90,90,-0.2,4.2]); 

xlabel('Input  Direction  of  Arrival  ( degrees )') 
grid 

%  Find  the  angle  from  the  3-bit  RNS  Modulus  code 
mslogl 

%  Change  3-bit  RNS  modulus  code  to  six-bit  binary  code 
bin_six 

%  Find  the  number  of  glitch 
glitch 

%  Plots  the  Resolved  DOA  and  Quantization  Error  as  a  Function  of  the  Input  Direction  of 

Arrival 

figure(2) 

whitebg 

plot(theta,th,'b') 

xlabel('Input  Direction  of  Arrival  ( degrees )') 
ylabel('Resolved  DOA  and  Quantization  Error  ( degrees  )') 
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title(['Encoding  Error  in  Quatization  (glitch=  ',num2str(g_count),', sampling 

distancerds- ,num2str(ds*3),'  degree)'])%,  ',num2str(l/ds),'samples  in  3  degree  )']) 

axis([-90  90  -9090]) 

grid 

print 

figure(3) 
whitebg 
subplot(21 1) 
plot(theta,modlp) 
axis([-90,90,-4,4]); 

xlabel('Input  Direction  of  Arrival  ( degrees )') 

title(['Mode  3  Foldingwave  and  Ou^ut  code  (sampling  distance:ds- ,num2str(ds*3),' 

degree)'])%,  ',num2str(l/ds), 'samples  in  3  degree )']) 

grid 

subplot(212) 

stairs(theta,modl) 

axis([-90,90,-0.2,2.2]); 

xlabel('Input  Direction  of  Arrival  ( degrees )') 
grid 

save  binl 

%  interpolate  the  Encoding  error 
interpl 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%°/»%%“/» 
%  This  program  'Treshold'  calculates  the  original  threshold  levels  for  setting  the  RNS  % 
%  folding  wave  circuit  in  the  0.3  degree  sampling  distance  % 


%Program  2 
%  thres.m 

global  ml 
global  m2 
global  m3 
global  fo 
global  jfr 
global  ds 
global  M 
global  L 

%  Calculates  the  Threshold  Values  for  Modulus  ml 
i=-l:2:l; 

Tl=(pi/3)*i; 

Tl=sprintfC%6.3f,Tl); 

Tl=sscanf(Ti;%f); 

%  Calculates  the  Threshold  Values  for  Modulus  m2 
i=-l:l; 

T2=(pi/2)*i; 

T2=sprintf('%6.3f,T2); 

T2=sscanf(T2;%f); 

%  Calculates  the  Threshold  Values  for  Modulus  m3 
i=-3:2:3; 

T3=(pi/5)*i; 

T3=sprintf('%6.3f,T3); 

T3=sscanf(T3,'%f); 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%This  program  'RNS_ANGLE'  is  to  determine  the  resolved  Direction  Of  Arrival  from  % 
%  the  RNS  output  in  the  moduli  set  (3,4,5)  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  program  3 
%  set  3-Bit  RNS  Code 
res=[modr  mod2'  mod3']; 

%  Number  of  Folds  Over  the  Entire  field  of  View  for  Modulus  ml 
iif=fr/fo*M/ml; 

%Total  Number  of  Quantization  Levels  Over  the  Entire  Field  of  View 
nql=nf'‘ml; 

%Generates  the  Quantization  Level  Transistion  Angles 
p=0:2/nql:l; 

p=[  -rot90(p(2:size(p,2)))'  p]; 
trang=l/rad*asin(p); 

%Calculates  the  Midvalue  of  Each  Quantization  Level 
for  n=l  :size(p,2)-l ; 

midval(n)=trang(n)+(trang(n+ 1  )-trang(n))/2 ; 
qlw(n)=midval(n)-trang(n); 
end 

if  rem(size(midval,2),2)=0; 
a=midval; 
b=size(a,2)/2; 
end 

if  rem(size(midval,2),2)=l; 
a=midval; 
b=(size(a,2)-i- 1  )/2; 
end 

if  size(a,2)“60; 
a=[zeros(l,b)  a  zeros(l,b)]; 
end 

b=size(a,2)/2; 

%Calculates  the  Resolved  Direction  of  Arrival 
foru=l:M/ds; 
if  res(u,:)=([0,0,0]); 
th(u)=a(l); 

elseif  res(u,:)=([  1,1,1]); 
th(u)=a(2); 
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elseif  res(u,:)==([2,2^]); 
th(u)=a(3); 

elseif  res(u,:)==([0,3,3]); 
th(u)=a(4); 

elseif  res(u,:)=([  1,0,4]); 
th(u)=a(5); 

elseif  res(u,:)=([2,l,0]); 
th(u)=a(6); 

elseif  res(u,:)=([0,2,l]); 
th(u)=a(7); 

elseif  res(u,:)==([l53,2]); 
th(u)=a(8); 

elseif  res(u,:)=([2,0,3]); 
th(u)=a(9); 

elseif  res(u,:)=([0,l,4]) ; 
th(u)=a(10); 

elseif  res(u,:)=(P  ,2,0]); 
th(u)=a(ll); 

elseif  res(u,:)=([2,3,l]); 
th(u)=a(12); 

elseif  res(ii,:)=([0,0,2]) ; 
th(u)=a(13); 

elseif  res(u,  :)==([  1,1, 3]); 
th(u)=a(14); 

elseif  res(u,:)=([2,2,4]) ; 
th(u)=a(15); 

elseif  res(u,;)=([0,3,0]); 
th(u)=a(16); 

elseif  res(u,;)==([l,0,l]); 
th(u)=a(17); 

elseif  res(u,;)=([2,l,2]) ; 
th(u)=a(18); 

elseif  res(u,:)=([0,2,3]); 
th(u)=a(19); 

elseif  res(u,:)=([l, 3,4]) ; 
th(u)=a(20); 

elseif  res(u,;)==([2,0,0]); 
th(u)=a(21); 

elseif  res(u,:)=([0,l,l]); 
th(u)=a(22); 

elseif  res(u,  :)=([!, 2,2]) ; 
th(u)=a(23); 


elseif  res(u,:)===([2,3,3]); 
th(u)=a(24); 

elseif  res(u,:)=([0,0,4]); 
th(u)=a(25); 

elseif  res(ii,:)==([l,l,0]); 
th(u)=a(26); 

elseif  res(xi,:)==([2,2,l]); 
th(u)=a(27); 

elseif  res(u,;)==([0,3,2]) ; 
th(u)=a(28); 

elseif  res(u,:)=([l,0,3]); 
th(u)=a(29); 

elseif  res(u,:)=([2, 1,4]) ; 
th(u)=a(30); 

elseif  res(u,:)=([0,2,0]); 
th(u)=a(31); 

elseif  res(u,  :)==([  1,3,1]); 
th(u)=a(32); 

elseif  res(u,:)— ([2,0,2]) ; 
th(u)=a(33); 

elseif  res(u,:)=([0,l,3]); 
th(u)=a(34); 

elseif  res(u,  :)==([  1,2,4])  ; 
th(u)=a(35); 

elseif  res(u,:)=([2,3,0]); 
th(u)=a(36); 

elseif  res(u,:)=([0,0,l]); 
th(u)=a(37); 

elseif  res(u,:)=([l,l, 2]) ; 
1h(u)=a(38); 

elseif  res(u,:)=([2,2,3]); 
th(u)=a(39); 

elseif  res(u,:)=([0,3,4]) ; 
th(u)=a(40); 

elseif  res(u,:)=([l,0,0]); 
th(u)=a(41); 

elseif  res(u,  :)=([2, 1,1]); 
th(u)=a(42); 

elseif  res(u,:)=([0,2,2]) ; 
th(u)=a(43); 

elseif  res(u,:)=([l,3,3]); 
th(u)=a(44); 


elseif  res(u,:)=([2,0,4]); 
th(u)=a(45); 

elseif  res(u,:)=([0,l,0]); 
th(u)=a(46); 

elseif  res(u,:)=([l,2,l]); 
th(u)=a(47); 

elseif  res(Uj:)=([2,3,2]) ; 
th(u)=a(48); 

elseif  res(u,:)=([0,0,3]); 
th(u)=a(49); 

elseif  res(u, :)=([!, 1,4]) ; 
th(u)=a(50); 

elseif  res(u,:)=([2,2,0]); 
th(u)=a(51); 

elseif  res(u,:)==([0,3,l]); 
th(u)=a(52); 

elseif  res(u, :)=([!, 0,2]) ; 
th(u)=a(53); 

elseif  res(u,:)=([2,l,3]); 
th(u)=a(54); 

elseif  res(u,:)==([0,2,4]) ; 
th(u)=a(55); 

elseif  res(u,:)=([l,3,0]); 
th(u)=a(56); 

elseif  res(u,:)=([2,0,l]); 
th(u)=a(57); 

elseif  res(u,:)=([0,l,2]) ; 
th(u)=a(58); 

elseif  res(u,:)=([l,2,3]); 
th(u)=a(59); 
else  res(u,:)=([2,3,4]) ; 
th(u)=a(60); 
end 
end 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"/«%%%% 
%  This  program  'SIX_BINARY'  converts  the  RNS  code  output  % 

%  to  6-bit  binary  code  at  the  instant  input  angle  % 

%  program  4 
%  bin_six.m 
%  make  six-binary  code 
for  u=l  :M/ds 
if  res(u,:)=([0,0,0]) 
thl(u)=l; 

binary(u,:)=([0  0  0  0  0  0]); 

elseif  res(u, :)=([!, 14]) 

thl(u)=2; 

binary(u,:)=([0  0  0  0  01]); 

elseif  res(u,:)— ([2,2,2]) 
thl(u)=3; 

binary(u,;)=([0  0  0  01  0]); 

elseif  res(u,:)=([0,3,3]) 
thl(u)=4; 

binary(u,:)=([0  0  0  01  1]); 

elseif  res(u,:)^'^^[  1  >0,4]) 
thl(u)=5; 

binary(u,:)=([0  0  0  1  0  0]); 

elseif  res(u,:)==([2,l,0]) 
thl(u)=6; 

binary(u,:)=([0  0  01  0  1]); 

elseif  res(u,;)=([0,2,l]) 
thl(u)=7; 

bmary(u,:)=([0  0  0  1  1  0]); 

elseif  res(u,:)==([l  ,3,2]) 
thl(u)=8; 

binary(u,:)=([0  0  0  1  1  1]); 
elseif  res(u,:)=([2,0,3]) 
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thl(u)=9; 

bmary(u,;)=([0  01  0  0  0]); 

elseif  res(u,:)=([0,l,4]) 
thl(u)=10; 

binary(u,:)=([0  01  0  0  1]); 

elseif  res(u,:) — ([1 ,2,0]) 
thl(u)=ll; 

binary(u,:)=([0  0101  0]); 

elseif  res(u,:)=([2,3,l]) 
thl(u)=12; 

binary(u,:)=([0  0101  1]); 
elseif  res(u,:)==([0,0,2]) 

thl(u)=13; 

binary(u,:)=([0  01  10  0]); 

elseif  res(u, :)=([!, 13]) 
thl(u)=14; 

binary(u,:)=([0  0110  1]); 

elseif  res(u,:)=([2,2,4]) 
thl(u)=15; 

binary(u,:)=([0  0  111  0]); 

elseif  res(u,:)=([0,3,0]) 
thl(u)=16; 

binary(u,:)=([0  0  111  1]); 

elseif  res(u, :)=([!, 0»1]) 
thl(u)=17; 

biiiary(u,:)=([0  1  0  0  0  0]); 

elseif  res(u,:)=([2,13]) 
thl(u)=18; 

binary(u,:)=([0  1  0  0  0  1]); 

elseif  res(u,:)=([0,23]) 
thl(u)=19; 

binary(u,:)=([0  10  01  0]); 


elseif  res(ii,  :)==([  1,3 ,4]) 
thl(u)=20; 

binary(u,:)=([0  10  01  1]); 

elseif  res(u,:)==([2,0,0]) 
thl(u)=21; 

binary(u,:)=([0  1010  0]); 

elseif  res(u,:)=([0,l,l]) 
thl(u)=22; 

binary(u,:)=([0  1010  1]); 

elseif  res(u,:)=([l,2,2]) 
thl(u)=23; 

bmaiy(u,:)=([0  1011  0]); 

elseif  res(u,:)=([2,3,3]) 
thl(u)=24; 

binary(u,:)=([0  10111]); 

elseif  res(u,:)=([0,0,4]) 
thl(u)=25; 

bmary(ii,:)=([0  1  1  0  0  0]); 

elseif  res(u,  :)=([  1  ?  1  jO]) 
thl(u)=26; 

binary(u,:)=([0  1  10  01]); 

elseif  res(u,:)=([2,2,l]) 
thl(u)=27; 

binary(u,:)=([0  1101  0]); 

elseif  res(u,:)==([0,3,2]) 
thl(u)=28; 

binary(u,:)=([0  110  11]); 

elseif  res(u,:)=([l,0,3]) 
thl(u)=29; 

binary(u,:)=([0  1110  0]); 


elseif  res(u,:)=([2,l,4]) 


thl(u)=30; 

binary(u,:)=([0  1  1  101]); 

elseif  res(u,:)=([0,2,0]) 
thl(u)=31; 

binary(u,:)=([0  1111  0]); 

elseif  res(u,:)=([l,34]) 
thl(u)=32; 

binary(u,:)=([0  1111  1]); 

elseif  res(u,:)=([2,0,2]) 
thl(u)=33; 

binary(u,:)=([l  0  0  0  0  0]); 

elseif  res(u,:)==([0,13]) 
thl(u)=34; 

binary(ii,  :)=([!  0  0  0  0  1]); 

elseif  res(u,;)=([l,2,4]) 
thl(u)=35; 

binaiy(u, :)=([!  0  0  01  0]); 

elseif  res(u,:)=([2,3,0]) 
thl(u)=36; 

binary(u,;)=([l  0  0  0  1  1]); 

elseif  res(u,:)=([0,0,l]) 
thl(u)=37; 

binary(u,:)=([l  0  01  0  0]); 

elseif  res(u, :)=([!  ,1 ,2]) 
thl(u)=38; 

binary(u,:)=([l  0  010  1]); 

elseif  res(u,:)=([2,2,3]) 
thl(u)=39; 

binary(u,:)=([l  0  011  0]); 

elseif  res(u,:)=([0,3,4]) 
thl(u)=40; 

binary(u, :)=([!  0  011  1]); 


elseif  res(u,;)=([l,0,0]) 
thl(u)=41; 

binary(u, :)=([!  010  0  0]); 

elseif  res(u,:)=([2,l,l]) 
thl(u)=42; 

binaiy(u,:)=([l  010  0  1]); 

elseif  res(u,:)==([0,2,2]) 
thl(u)=43; 

binary(u, :)=([!  0101  0]); 

elseif  res(xi,:)==([l,353]) 
thl(u)=44; 

binaiy(u,:)=([l  0101  1]); 

elseif  res(u,:)=([2,0,4]) 
thl(u)=45; 

bmary(u,:)=([l  0110  0]); 

elseif  res(u,:)=([0,l,0]) 
thl(u)=46; 

binary(u, :)=([!  0110  1]); 

elseif  res(u, ;)=([!, 2,1]) 
thl(u)=47; 

binary(u, :)=([!  0111  0]); 

elseif  res(u,:)=([2,3,2]) 
thl(u)=48; 

binaiy(ii,;)=([l  0111  1]); 

elseif  res(u,:)=([0,0,3]) 
thl(u)=49; 

binaiy(u,  :)=([!  1  0  0  0  0]); 

elseif  res(u, :)==([!, 1,4]) 
thl(u)=50; 

binary(u,:)=([l  10  0  0  1]); 


elseif  res(vi,:)=([2,2,0]) 


thl(u)=51; 

binary(u,;)=([l  10  01  0]); 

elseif  res(u,;)==([0,3,l]) 
thl(u)=52; 

binary(i4:)=([l  10  01  1]); 

elseif  res(u,:)=([l,0,2]) 
thl(u)=53; 

binary(u,:)=([l  1010  0]); 

elseif  res(u,:)=([2,l,3]) 
thl(u)=54; 

binary(u, :)=([!  1010  1]); 

elseif  res(u,:)=([0,2,4]) 
thl(u)=55; 

binary(u,:)=([l  1011  0]); 

elseif  res(u,:)=([l,3,0]) 
thl(u)=56; 

binary(u, :)=([!  1011  1]); 

elseif  res(u,:)=([2,0,l]) 
thl(u)=57; 

binary(u,:)=([l  1  1  0  0  0]); 

elseif  res(u,:)=([0,l,2]) 
thl(u)=58; 

binary(u, :)=([!  110  0  1]); 

elseif  res(u,:)— ([1,2,3]) 
thl(u)=59; 

binary(u,;)=([l  110  1  0]); 

elseif  res(u,:)=([2,3,4]) 
thl(u)=60; 

binary(u,:)=([l  1101  1]); 

end 

end 

%binary' 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  This  program  'Shift_LSB'  performs  the  interpolation  using  Shift  LSB  Method  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%“/«% 

%  program  5 
%  inteiplm 

%  set  the  six-addittional  threshol 
adth3 

%  This  is  the  operating  loop  until  'M/ds'  number  of  parity  samples 
i=l; 

while  i<M/ds 

a=l;  %  initial  value  of 'a'  is  1 

%  number  of  continuous  error  parity  in  5  parity  samples 
%  for  example,  when  error  is  1,  a=2 
while  parity(i)~=l  %  to  find  the  first  T  parity(error  parity) 
i=i+l ;  %  i  is  the  location  of  error 

end 

while  parity(i)-^0  %  to  find  the  next  error  parity('O'value) 
i=i+l; 
end 

while  parity(i)~^l  %  to  find  the  continuous  number  of  error  parity 
i=i+l; 
a=a+l; 

if  i>=M/ds+l  ,break,end  %  if  the  number  of  parity  is 
%  bigger  than  M/ds, then  finish  loop  ! 
end 

if  i>=M/ds+ 1  ,break,end 
i=i-a+l;  %  real  location  value  of 'O' parity 
b=l;  %  highest  level  of  LSB  value 
while  binary(i-l,b)=binary(i+a-l,b) 
b=b+l; 
ifb=7 
break 
end 
end 
ifb=l 
for  k=i+3:-l:i 
forj=l:6 

binary(ko)=binary(k-l  ,j); 
end 
end 
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end 

%  when  the  number  of  parity  error  is  one  (a=2) 
ifa=2 

%  change  the  MSB  value 
forj=l:b-l 

binary(i+3j)=binary(i+2j); 
binaiy(i+2  j)=binary(i+l  J); 
binary(i+l  J)=binary(i+1  J); 
binary(  i,  j)=binary(i-l  j); 
end 

%  change  the  LSB  value 
if  b<7 

fork=i+3:-l:i 

forj=b;6 

binary(k  j)=binary(k-l  j); 
end 
end 
end 

%  when  the  number  of  parity  error  is  two  (a=3) 
elseif  a==3 

%  change  the  MSB  value 
forj=l:b-l 

binary(i+3  j)=binary(i+2  j); 
binary(i+2  j)=binary(i-l  j); 
binary(i+lj)=binary(i-lj); 
binary(  i,  j)=binary(i-l  j); 

end 

%  change  the  LSB  value 
if  b<7 

for  k=i+3:-l:i 
forj=b:6 

binary(k  j)=binary(k-l  j); 
end 
end 
end 

%  when  the  number  of  parity  error  is  three  (a=4) 
elseif  a=4 

%  change  the  MSB  value 
for  j=l:b-l 
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binary(i+3j)=binary(i-lj); 
binary0+2j)=binaiy(i-lj); 
binary(i+l  j)=binaiy(i-l  j); 
binary(  i,  j)=binary(i-l  j); 
end 

%  change  the  LSB  value 
ifb<7 

fork=i+3:-l:i 

forj=b:6 

binary  (kj)=binary(k- 1  j); 
end 
end 
end 

%  when  the  number  of  parity  error  is  four 
%a=5 

%  change  the  all  binary  value 
else 

for  k?=i+3:-l  :i+l 
forj=l:6 

binary(kj)=binary(k- 1  J); 
end 
end 

end  %  if  a==2 

%  printout  6-bit  binary-code  after  interpolation 
if(i+3)<(M/ds) 
v=i-l:i+3; 
parity(v) 
binary  (v,;)' 
end 
i=i+a; 

end  %  while  i<=M/ds 

%  change  6-bit  binary  code  to  Modulus  (ml,ni2,m3) 
ant_bin 

%  Calculates  the  Resolved  Direction  of  Arrival  after  interpolation 
nislog2 

%  find  the  number  of  glitch 
gi_count=0; 
for  u=l  :M/ds-l 
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max(u+l  )=th2(u+l  )-th2(u); 
if  u<M/(2*ds) 
if  max(u+l)>  1 1 

parity(u+l  )=0;  %  possible  bad  samples 

gi_count=gi_coxmt+l ; 
u+l4nax(u+l) 
parity(u:u+4) 
binary(u:u+4,:)' 
else 

parity(u+l)=l ;  %  sample  is  good 
end 
else 

if  max(u+l)>  1 1 

parity(u)=0;  %  possible  bad  samples 
gi_coimt=gi_count+l ; 
u+l,max(u+l) 
parity(u-r.u+3) 
binary(u-l:u+3,:)' 
else 

parity(u)=l ;  %  sample  is  good 
end 

end%  ifu<M/(2*ds) 
end  %  for 


figure(5) 

whitebg 

plot(theta,th2,'b') 

xlabel('Input  Direction  of  Arrival  ( degrees )') 

ylabel('Resolved  DOA  and  Quantization  Error  ( degrees )') 

title(['Interpolation  Result,' ,num2str(LSB_percent),'  %  LSB-Shift  Method  (  glitch= 

',num2str(gi_count),',  sampling  distancerds- ,num2str(ds*3),'  degree)']) 

axis([-90  90  -90  90]) 

save  bin2  Tl_add  gi_cormt  binary  parity; 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  This  program  'ADD_THRESH'  calculates  the  additional  threshold  levels  % 

%  to  remove  the  quantization  error  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  program  6 
%  adth3.m 

global  ml 
global  m2 
global  m3 
global  fo 
global  fr 
global  ds 
global  rad 
global  M 
global  theta 
global  thet 

%LSBjpercent  =input('What  is  the  percentage  of  LSB  ?(%)') 
incre=0.6*LSB_percenV100;  %  Set  the  threshold  by  1  degree 

%  set  the  additional  thresholds 
Tl_add(  1  )=2*pi*incre/9-pi; 

Tl_add(2)=2*pi*(3-incre)/9-pi; 

Tl_add(3)=2*pi*(3+incre)/9-pi; 

Tl_add(4)=2*pi*(6-incre)/9-pi; 

Tl_add(5)=2*pi*(6+incre)/9-pi; 

Tl_add(6)=2*pi*(9-incre)/9-pi; 

Tl_add 

%  calculate  the  six-comparator  output 
for  i=l  :M/ds 

if  modlp(i)  <=  Tl_add(l) 
count(i)  =  0;  %  Bad  sample 

elseif  modlp(i)  <=  Tl_add(2) 
count(i)  =  1 ;  %  Good  sample 

elseif  modlp(i)  <=  Tl_add(3) 
count(i)  =  2; 

elseif  modlp(i)  <=  Tl_add(4) 
count(i)  =  3; 

elseif  modlp(i)  <=  Tl_add(5) 
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coimt(i)  =  4; 

elseif  modlp(i)  <=  Tl_add(6) 
comit(i)  =  5; 
else 

count(i)  =  6; 
end 
end 


%  parity  cicuit 
for  i=l:M/ds 
if  rem(count(i),2)  =  0 
parity(i)=0;  %  possible  bad 
else 

parity(i)=  1 ;  %  good  parity 

end 
end 

figure(4) 
whitebg 
subplot(21 1) 
plot(theta,modlp) 
axis([-90,90,-4,4]); 

xlabel('Input  Direction  of  Arrival  ( degrees )') 
grid 

subplot(212) 
stairs(thet, count) 
axis([-90,90r0.2,6.2]); 

xlabel('Input  Direction  of  Arrival  (  degrees  )') 

grid 

print 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  This  program  'ANTJBINARY*  converts  the  six-bit  binary  code  % 

%  to  the  RNS  code  after  interpolation  the  encoding  error  % 

%%o/o%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  program  7 
%  ant_bin.m 

%  change  six-binary  code  to  RNS  Code 
foru=l:M/ds 

if  binary(u,:)=([0  0  0  0  0  0]) 
res2(u,:)=([0,0,0]); 
thl(u)=l; 

elseif  binary(u,  :)=([0,0,0,0,0, 1  ]) 
res2(u,:)=([l,l,l]); 

thl(u)=2; 

%binary(u,:) 

%res(u,:) 

elseif  binary(u,:)=([0  0  0  010]) 
res2(u,:)=([2,2,2]); 
thl(u)=3; 

elseif  binaiy(u,:)==([0  0  0  01  1]) 
res2(u,:)=([0,3,3]); 
thl(u)=4; 

elseif  binary(u,:)=([0  0  01  0  0]) 

res2(u,:)=([U0,4]); 

thl(u)=5; 

elseif  binaiy(u,:)==([0  0  0101]) 
res2(u,:)=([2,l,0]); 
thl(u)=6; 

elseif  binary(u,:)=([0  0  01  1  0]) 
res2(u,:)=([0,2,l]); 
thl(u)=7; 

elseif  binary(u,:)=([0  0  01  1  1]); 
res2(u,:)=([l,3,2]); 
thl(u)=8; 

elseif  binary(u,;)=([0  01  0  0  0]); 
res2(u,:)=([2,0,3]); 
thl(u)=9; 

elseif  binary(u,:)=([0  01  0  0  1]); 
res2(u,:)=([0,l,4]) ; 
thl(u)=10; 

elseif  binary(u,:)=([0  0  10  1  0]); 
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res2(u,:)=([l,2,0]); 

thl(u)=ll; 

elseif  binary(u,:)=([0  0101  1]); 

res2(u,:)=([23,l]); 

thl(u)=12; 

elseif  binary(u,:)==([0  01  10  0]); 
res2(u,:)=([0»0,2]) ; 
thl(u)=13; 

elseif  binary(u,:)=([0  01  101]); 

res2(u,:)=([l,13]); 

thl(u)=14; 

elseif  binary(u,:)=([0  0111  0]); 
res2(u,:)=([2,2,4]) ; 
thl(u)=15; 

elseif  binary(u,:)==([0  0111  1]); 
res2(u,:)=([0,3,0]); 
thl(u)=16; 

elseif  binary(u,:)=([0  1  0  0  0  0]); 

res2(u,:>=([lAl]); 

thl(u)=17; 

elseif  binary(u,:)=([0  1  0  0  0  1]); 
res2(u,:)=([2,l,2]) ; 
thl(u)=18; 

elseif  binary(u,:)==([0  10  01  0]); 
res2(u,:)=([0,2,3]); 
thl(u)=19; 

elseif  bmary(u,:)==([0  10  01  1]); 
res2(u,:)=([l,3,4]) ; 
thl(u)=20; 

elseif  binary(u,:)==([0  1010  0]); 

res2(u,:)=([2,0,0]); 

thl(u)=21; 

elseif  bmary(u,:)=([0  10101]); 
res2(u,:)=([0,l,l]); 
thl(u)=22; 

elseif  binary(u,:)=([0  1011  0]); 
res2(u,:)=([l,2,2]) ; 
thl(u)=23; 

elseif  binary(u,:)==([0  1011  1]); 
res2(u,:)=([2,3,3]); 
thl(u)=24; 

elseif  binary(u,:)=([0  1  1  0  0  0]); 


res2(u,:)=([0,0,4]); 

thl(u)=25; 

elseif  binary(u,;)=([0  110  0  1]); 
res2(u,:)=([14,0]); 
thl(u)=26; 

elseif  binary(u,:)=([0  1101  0]); 
res2(u,:)=([2,2,l]); 
thl(u)=27; 

elseif  binaiy(i4:)=([0  1101  1]); 
res2(u,:)=([0,3,2]) ; 
thl(u)=28; 

elseif  binary(u,:)=([0  1110  0]); 
res2(u,:)=([l,0,3]); 
thl(u)=29; 

elseif  binary(u,:)=([0  1110  1]); 
res2(u,:)=([2,l,4]) ; 
thl(u)=30; 

elseif  bmary(u,:)==([0  1111  0]); 
res2(u,:)=([0,2,0]); 
thl(u)=31; 

elseif  binary(u,:)==([0  1111  1]); 

res2(u,:)=([l,34]); 

thl(u)=32; 

elseif  binary(u,:)=([l  0  0  0  0  0]); 
res2(u,:)=([2,0,2]) ; 
thl(u)=33; 

elseif  bmaiy(u,:)=([l  0  0  0  0  1]); 

res2(u,:)=([0,l,3]); 

thl(u)=34; 

elseif  biiiary(u,:)=([l  0  0  0  1  0]); 
res2(u, :)=([!  »2, 4])  ; 
thl(u)=35; 

elseif  binary(ii,  :)==([ 1  0  0  0  1  1]); 
res2(u,:)=([2,3,0]); 
thl(u)=36; 

elseif  binary(u,:)=([l  0  0  1  0  0]); 

res2(u,:)=([0,04]); 

thl(u)=37; 

elseif  binaiy(u,:)=([l  0  010  1]); 
res2(u,:)=([l,l,2]) ; 
thl(u)=38; 

elseif  binary(u,:)=([l  0  011  0]); 


res2(u,:)=([2,2,3]); 

thl(u)=39; 

elseif  binary(u,:)=([l  0011  1]); 
res2(u,:)=([03,4]) ; 
thl(u)=40; 

elseif  binary(u,:)—([l  01  0  0  0]); 
res2(u,:)=([lA0]); 
thl(u)=41; 

elseif  binary(u,:)==([l  010  0  1]); 
res2(u,:)=([2,l,l]); 
thl(u)=42; 

elseif  bmary(u,  :)=([!  0101  0]); 
res2(u,:)=([0,2^]) ; 
thl(u)=43; 

elseif  binary(u, :)=([!  0101  1]); 
res2(u,:)=([13,3]); 
thl(u)=44; 

elseif  bmary(u, :)==([!  0110  0]); 
res2(u,:)=([2,0,4]); 
thl(u)=45; 

elseif  binary(u, :)=([!  01  101]); 
res2(u,;)=([0,l,0]); 
thl(u)=46; 

elseif  binary(u,  ;)==([  1  0111  0]); 
res2(u, :)=([!  ,2,1]); 
thl(u)=47; 

elseif  binaiy(u,:)==([l  0111  1]); 
res2(u,:)=([2,3,2]) ; 
thl(u)=48; 

elseif  binary(u,  :)=([!  1  0  0  0  0]); 
res2(u,:)=([0,0,3]); 
thl(u)=49; 

elseif  binaiy(u,:)=([l  10  0  0  1]); 
res2(u,:)=([14,4]) ; 
thl(u)=50; 

elseif  binary(u,:)=([l  10  01  0]); 
res2(u,:)=([2,2,0]); 
thl(u)=51; 

elseif  binary(u,;)=([l  10  01  1]); 
res2(u,:)=([0,3,l]); 
thl(u)=52; 

elseif  binary(u,:)=([l  1010  0]); 


res2(u,:)=([l,0,2]) ; 
thl(u)=53; 

elseif  bmary(u,:)=([l  1010  1]); 

res2(u,:)=([2,l,3]); 

thl(u)=54; 

elseif  binary(u,:)=([l  1011  0]); 
res2(u,:)=([0,2,4]) ; 
thl(u)=55; 

elseif  binaiy(u,  :)==([  1  1011  1]); 
res2(u,:)=([13,0]); 
thl(u)=56; 

elseif  binary(u,:)=([l  1  1  0  0  0]); 
res2(u,:)=([2,0,l]); 
thl(u)=57; 

elseif  binary(u,:)=([l  110  0  1]); 
res2(ii,:)=([0,l,2]) ; 
thl(u)=58; 

elseif  binary(u,  :)=([!  1101  0]); 

res2(u,:)=([l,2,3]); 

thl(u)=59; 

elseif  binary(u,  :)==([  1  110  1  1]); 

res2(u,:)=([2,3,4]) ; 

thl(u)=60; 

end 

end 


%binary' 


%Vo%%%Vo%VoVoVoVoVoVo%Vo%%%%%%%%%%%%%Vo%%%VtP/o%%VoVo%%%% 

%  This  program  'RNS_angle2'  is  to  determine  the  resolved  Direction  Of  Arrival  % 
%  the  RNS  output  in  the  moduli  set  (3,4,5)  % 


%  Program  8 
%  mslog2.m 

%  Number  of  Folds  Over  the  Entire  field  of  View  for  Modulus  ml 
nf=fr/fo*M/ml; 

%Total  Number  of  Quantization  Levels  Over  the  Entire  Field  of  View 
nql=nf*ml; 

%Generates  the  Quantization  Level  Transistion  Angles 
p=0:2/nql:l; 

p=[  -rot90(p(2:size(p,2)))'  p]; 
trang=l/rad*asin(p); 

%Calculates  the  Midvalue  of  Each  Quantization  Level 
forn=l:size(p,2)-l 

midval(n)=trang(n)+(trang(n+l)-trang(n))/2; 

%qlw(n)=midval(n)-trang(n); 

end 

if  rem(size(midval,2),2)=0 

a=midval; 

b=size(a,2)/2; 

end 

if  rem(size(midval,2),2)==l 

a=midval; 

b=(size(a,2)+l)/2; 

end 

if  size(a,2)~=60 
a=[zeros(l,b)  a  zeros(l,b)]; 
end 

b=size(a,2)/2; 

%Calculates  the  Resolved  Direction  of  Arrival 
for  u=l:M/ds 
if  res2(u,;)=([0,0,0]); 
th2(u)=a(l); 

elseif  res2(u,:)=([l,l,l]); 
th2(u)=a(2); 

elseif  res2(u,:)=([2,2,2]); 
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th2(u)=a(3); 

elseif  res2(u,;)==([0,3,3]); 
th2(u)=a(4); 

elseif  res2(u,:)=([l,0,4]); 
th2(u)=a(5); 

elseif  res2(u,:)=([2,l,0]); 
th2(u)=a(6); 

elseif  res2(u,:)=([0,2,l]); 
th2(u)=a(7); 

elseif  res2(u,:)=([l,3»2]); 
th2(u)=a(8); 

elseif  res2(u,:)=([2,0,3]); 
th2(u)=a(9); 

elseif  res2(u,:)=([0,l»4]) ; 

th2(u)=a(10); 

elseif  res2(u,:)=([l^,0]); 

th2(u)=a(l  1); 

elseif  res2(u,:)==([2,3,l]); 

th2(u)=a(12); 

elseif  res2(u,:)==([0,0^]) ; 

th2(u)=a(13); 

elseif  res2(u,:)=([l,l,3]); 

th2(u)=a(14); 

elseif  res2(u,:)=([2,2,4]) ; 

th2(u)=a(15); 

elseif  res2(u,  :)=([0,3,0]); 

th2(u)=a(16); 

elseif  res2(u,:)=([l,0,l]); 

tli2(u)=a(17); 

elseif  res2(u,:)=([2, 1,2]) ; 

th2(u)=a(18); 

elseif  res2(u,:)=([0,2,3]); 

th2(u)=a(19); 

elseif  res2(u,:)=([l,3,4]) ; 

th2(u)=a(20); 

elseif  res2(u,:)=([2,0,0]); 

th2(u)=a(21); 

elseif  res2(u,:)=([0,l,l]); 

th2(u)=a(22); 

elseif  res2(u,:)=([  1,2, 2]) ; 

th2(u)=a(23); 

elseif  res2(u,:)=([2,3,3]); 


th2(u)=a(24); 

elseif  res2(ii,:)=([0,0,4]); 

th2(u)=a(25); 

elseif  res2(u,:)==([l,l,0]); 

th2(u)=a(26); 

elseif  res2(u,:)=([2,2,l]); 

th2(u)=a(27); 

elseif  res2(u,:)=([0,3,2]) ; 

th2(u)=a(28); 

elseif  res2(u,:)=([l,0,3]); 

th2(u)=a(29); 

elseif  res2(u,:)=([2, 1,4]) ; 

th2(u)=a(30); 

elseif  res2(u,:)==([0,2,0]); 

th2(u)=a(31); 

elseif  res2(u,:)==([l,3,l]); 

th2(u)=a(32); 

elseif  res2(u,:)=([2,0,2]) ; 

th2(u)=a(33); 

elseif  res2(u,:)=([0,l,3]); 

th2(u)=a(34); 

elseifres2(u,:)===([l,2,4]) ; 

th2(u)=a(35); 

elseif  res2(u,:)==([2,3,0]); 

th2(u)=a(36); 

elseif  res2(u,:)=([0,0,l]); 

th2(u)=a(37); 

elseif  res2(u,:)==([l,  1,2]) ; 

1h2(u)=a(38); 

elseif  res2(u,:)=([2,2,3]); 

th2(u)=a(39); 

elseif  res2(u,:)==([0,3,4]) ; 

th2(u)=a(40); 

elseif  res2(u,:)=([l,0,0]); 

th2(u)=a(41); 

elseif  res2(u,:)=([2,l,l]); 

th2(u)=a(42); 

elseif  res2(u,:)==([0,2,2]) ; 

th2(u)=a(43); 

elseif  res2(u,:)==([l,3,3]); 

th2(u)=a(44); 

elseif  res2(u,:)=([2,0,4]); 


th2(u)=a(45); 
elseif  res2(u,:)=([0,l,0]); 
th2(u)^a(46)  j 
elseif  res2(u,:)=([l,2,l]); 
th2(u)=a(47); 
elseif  res2(u,:)=([2,3,2]) ; 
th2(u)=a(48); 
elseif  res2(u,;)=([0,0,3]); 
th2(u)=a(49); 
elseif  res2(u,:) — ([1,1,4]) ; 
th2(u)=a(50); 
elseif  res2(u,:)=([2,2,0]); 
th2(u)=a(51); 
elseif  res2(u,:)=([0,3,l]); 
11i2(u)=a(52); 
elseif  res2(u,:)=([l,0, 2]) ; 
th2(u)=a(53); 
elseif  res2(u,:)=([2,l,3]); 
th2(u)=a(54); 
elseif  res2(u,:)=([0, 2,4]) ; 
th2(u)=a(55); 
elseif  res2(u,:)==([l,3,0]); 
th2(u)=a(56); 
elseif  res2(u,:)=([2,0,l]); 
th2(u)=a(57); 
elseif  res2(u,:)=([0, 1,2]) ; 
th2(u)=a(58); 
elseif  res2(u,  :)=([!, 2,3]); 
th2(u)=a(59); 
else  res2(i4;)=([2,3,4]) ; 
th2(u)=a(60); 
end 
end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  This  program  'Last_Good'  performs  the  interpolation  using  Last  good-sample  shift  % 

%  program  9 
%  int_copy 

%  set  the  six-addittional  threshol 
adthS 

%  This  is  the  operating  loop  until  'M/ds'  number  of  parity  samples 

i=i; 

while  i<M/ds 

a=l;  %  initial  value  of 'a'  is  1 

%  number  of  continuous  error  parity  in  5  parity  samples 
%  for  example,  when  error  is  1,  a=2 
while  parity(i)~=l  %  to  find  the  first 'T  parity(error  parity) 
i=i+l ;  %  i  is  the  location  of  error 

end 

while  parity(i)~=0  %  to  find  the  next  error  parity('0'value) 
i=i+l; 
end 

while  parity(i)~^l  %  to  find  the  continuous  number  of  error  parity 
i=i+l; 
a=a+l; 

if  i>=M/ds+ 1 , break, end  %  if  the  number  of  parity  is 
%  bigger  than  M/ds,then  finish  loop  ! 
end 

if  i>=M/ds+l  ,break,end 

i=i-a+l;  %  real  location  value  of 'O' parity 

b=  1 ;  %  highest  level  of  LSB  value 

%  Copy  the  last  good  sample  to  error  samples 
for  k=i+a:-l  :i 
forj=l:6 

binary(k,j)=binary(i-2J); 
end 
end 
i=i+l ; 

end  %  while  i<=M/ds 

%  change  6-bit  binary  code  to  Modulus  (ml, m2, m3) 
ant  bin 
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%  Calcvilates  the  Resolved  Direction  of  Arrival  after  interpolation 
mslog2 

%  find  the  number  of  glitch 
%parit5r=ones(  1  ,M/ds); 
gi_co;mt=0; 
for  u=l  :M/ds-l 
max(u+ 1  )=th2(u+ 1  )-th2(u) ; 
if  u<M/(2*ds) 
if  max(u+l)>  1 1 

parity(u+l  )=0;  %  possible  bad  samples 

gi_count=gi_count+l ; 
u+l,max(u+l) 
parity(u:u+4) 
binary(u:u+4,:)' 
else 

parity(u+l  )=1 ;  %  sample  is  good 
end 
else 

ifmax(u+l)>  11 

parity(u)=0;  %  possible  bad  samples 
gi_count=gi_count+l ; 
u+l,max(u+l) 
parity(u-l  :u+3) 
binary(u-l:u+3,:)' 
else 

parity(u)=  1 ;  %  sample  is  good 
end 

end%  ifu<M/(2*ds) 
end  %  for 

figure(5) 

whitebg 

plot(theta,th2,'b') 

xlabel('Input  Direction  of  Arrival  ( degrees  )’) 
ylabel('Resolved  DOA  and  Quantization  Error  (  degrees  )') 

title(['Interpolation  Result,', num2str(LSB_percent);  %  Shift  Last-Good  Bit  Method  (gntch= 
%',num2str(gi_count),',  sampling  distance:ds- ,num2str(ds*3),'  degree)']) 
axis([-90  90  -90  90]) 
print 

save  bin3  Tl_add  gi_count  binary  parity; 


%  This  program  'RAND_LSB'  performs  the  interpolation  using  random  LSB  Method 
% 

%program  10 
%  int_rand 

%  set  the  six-addittional  threshol 
adth3 

%  This  is  the  operating  loop  until  'M/ds'  number  of  parity  samples 

i=i; 

while  i<M/ds 

a=l;  %  initial  value  of 'a'  is  1 

%  number  of  continuous  error  parity  in  5  parity  samples 
%  for  example,  when  error  is  1,  a=2 

while  parity(i)~=l  %  to  find  the  first 'T  parity(error  parity) 
i=i+l ;  %  i  is  the  location  of  error 

end 

while  parity(i)~^0  %  to  find  the  next  error  parity('0'value) 
i=i+l; 
end 

while  parity(i)~=l  %  to  find  the  continuous  number  of  error  parity 
i=i+l; 
a=a+l; 

if  i>=M/ds+l , break, end  %  if  the  number  of  parity  is 
%  bigger  than  M/ds,then  finish  loop  ! 

end 

if  i>=M/ds+ 1  ,break,end 
i=i-a+l;  %  real  location  value  of 'O' parity 
b=l ;  %  highest  level  of  LSB  value 
while  bmary(i-l,b)==binary(i+a-l,b) 
b=b+l; 
ifb=7 
break 
end 
end 
ifb=l 
fork=i+3;-l;i 

binary(k,:)=rem(randperm(6),2); 

end 
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end 


%  when  the  number  of  parity  error  is  one  (a=2) 
ifa=2 

%  change  the  MSB  value 
forj=l:b-l 

binary(i+3j)=binary(i+2j); 
binary(i+2  J)=binary(i+1  J); 
binary(i+l  j)=binary(i+l  J); 
binary(  i,  j)=binary(i- 1  j); 
end 

%  randomize  the  LSB  value 
if  b<7 

fork=i+3:-l;i 

binary(k,  :)=Tem(randperm(6),2); 
end 
end 

%  when  the  number  of  parity  error  is  two  (a=3) 
elseif  a=3 

%  change  the  MSB  value 
forj=l:b-l 

binary(i+3j)=binary(i+2j); 
binary(i+2J)=binary(i-lj); 
binary(i+l  J)=binary(i- 1  J); 
binary(  i,  j)=binary(i-l  j); 
end 

%  change  the  LSB  value 
if  b<7 

for  k=i+3;-l  :i 
forj=b:6 

binary(k  j)=binary(k- 1  J); 
end 
end 
end 

%  when  the  number  of  parity  error  is  three  (a=4) 
elseif  a==4 

%  change  the  MSB  value 
forj=l:b-l 

binary(i+3  J)=binary(i- 1  J); 
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binary(i+2  j)=binaiy(i-l  j); 
binary(i+l  j)=binary(i-l  J); 
binary(  i,  j>=binary(i-l  J); 
end 

%  change  the  LSB  value 
ifb<7 

fork=i+3:-l:i 

binary(k,:)=rem(randpemi(6),2); 

end 

end 

%  when  the  number  of  parity  error  is  four 
%  a=5 

%  change  the  all  binary  value 
else 

fork=i+3:-l:i+l 

forj=l:6 

binary(ky)=binary(k-lj); 

end 

end 

end  %ifa=2 

%  printout  6-bit  binary-code  after  interpolation 
i 

if(i+3)<(M/ds) 
v=i-l  :i+3; 
parity(v) 
binary(v,:)' 
end 
i=i+a; 

end  %  while  i<=M/ds 

%  change  6-bit  binary  code  to  Modulus  (ml, m2, m3) 
antbin 

%  Calculates  the  Resolved  Direction  of  Arrival  after  interpolation 
mslog2 

%  find  the  number  of  glitch 
parity=ones(  1  ,M/ds); 
gi_count=0; 
for  u=l:M/ds-l 
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max(u+ 1  )=th2(u+ 1  )-th2(u) ; 
ifu<M/(2*ds) 
ifmax(u+l)>  11 

parity(u+l)=0;  %  possible  bad  samples 
gi_count=gi_count+l ; 
u+l,max(u+l) 
parity(u:u+4) 
binary(u:u+4,:)' 
else 

parity(u+l)=l ;  %  sample  is  good 
end 
else 

ifmax(u+l)>  11 

parity(u)=0;  %  possible  bad  samples 
gi_count=gi_count+l ; 
u+l,max(u+l) 
parity(u-l:u+3) 
binary  (u-l:u+3,:)' 
else 

parity(u)=l ;  %  sample  is  good 

end 

end%  if u<M/(2*ds) 
end  %  for 

figure(5) 

whitebg 

plot(theta,th2,'b') 

xlabel('Input  Direction  of  Arrival  ( degrees  )') 
ylabel('Resolved  DOA  and  Quantization  Error  ( degrees )') 

title(['Interpolation  Result, ',num2str(LSB_percent),'  %  Random-LSB  Method  ( glitch= 
%',num2str(gi_count),’,  sampling  distancerds- ,num2str(ds*3),'  degree)']) 
axis([-90  90  -90  90]) 

%print 

save  bin4  Tl_add  gi_coimt  binary  parity; 


98 


LIST  OF  REFERENCES 


1.  Unnikrishna,  Pillai,  S.,  “A  New  Approach  to  Array  Geometry  for  Improved  Spatial 
Spectrum  Estimation,”  IEEE,  vol.  73,  no.  10,  October  1985. 

2.  Rodrigues,  Louis  and  E.  Moita,  “High-resolution  Residue  Antenna  Architectures  for 
Wideband  Direction  Finding,”  Master  Thesis,  Naval  Postgraduate  School,  June 
1996. 

3.  Pace,  P.  E.,  L.  E.  M.  Rodrigues  and  D.  C.  Jenn,  “High  Resolution  Residue  Antenna 
Architectures  for  Wideband  Direction  Finding,  Department  of  Electrical  and 
Computer  Engineering,”  submitted  for  presentation  at  the  1997  National  Radar 
Conference,  June  1996. 

4.  Travers,  D.  N.,  and  S.  M.  Hixon,  Abstracts  of  Available  Literature  on  Radio 
Direction  Finding,  Southwest  Research  Institute,  1966. 

5.  Benoit,  R.  C.  Jr.,  and  F.  Coughlin  Jr.,  “Designing  RDF  Antennas,”  Electronic 
Industries,  vo\.\Z,  pp.77-83,  April  1959. 

6.  Jenkins,  Hemodon  H.,  Small  Aperture  Radio  Direction  Finding,  Norwood,  MA: 
Artech  House  Inc.,  1991. 

7.  Simon,  Haykin,  Radio  Array  Signal  Processing  for  Angle  of  Arrival  Estimation, 
Prentice-Hdl,  Inc.,  Englewood  Cliffs,  New  Jersey,  1985. 

8.  Grosswald,  E.,  Topics  From  the  Theory  of  Numbers,  McMillan,  New  York,  1966. 

9.  Gething,  P.  J.  D.,  “High-Frequency  Direction  Finding,”  Proceedings  of  the  IEEE, 
vol.  113,  pp.  49-61,  January  1966. 

10.  Rhodes,  D-P..,  Introduction  to  Monopulse  Radar,  McGraw-Hill,  1959. 

11.  Frank,  Jay,  ed.,  IEEE  Standard  Dictionary  of  Electrical  and  Electronic  Terms 
ANSFIEEE  Standard,  Institute  of  Electrical  Engineers,  Inc.,  New  York,  1988. 


99 


100 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 

1 .  Defense  Technical  Information  Center  2 

8725  John  J.  Kingman  Rd.,  STE  0944 

Ft.  Belvoir,  VA  22060-6218 

2.  Dudley  Knox  Library  2 

Naval  Postgraduate  School 

41 1  Dyer  Rd. 

Monterey,  California  93943-5101 

3.  Chairman,  Code  EW  1 

Department  of  Electronic  Warfare  Academic  Group 

Naval  Postgraduate  School 
Monterey,  California  93943-5121 

4.  Professor  David  C.  Jeim,  Code  EC/Jn  1 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5121 

5.  Professor  Phillip.  E.  Pace,  Code  EC/Pc  3 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5 121 

6.  Space  and  Naval  Warfare  System  Command  1 

Department  of  the  Navy 

PMW-163 

Attn:  CAPT.  Coimell 
Washington,  DC  20363-5100 

7.  Space  and  Naval  Warfare  System  Command  1 

Department  of  the  Navy 

PMW-163 

Attn:  Mike  Wilson 

Washington,  DC  20363-5 1 00 


101 


8.  Space  and  Naval  Warfare  System  Command 
Department  of  the  Navy 

PMW-163 

Attn:  Ralph  Skiano 

Washington,  DC  20363-5100 

9.  Space  and  Naval  Warfare  System  Command 
Department  of  the  Navy 

PMW-163 

Attn:  CAPT.  Ristorcelli 
Washington,  DC  20363-5100 

10.  Byeong-Jun  Park 
368-14  Sansu-Dong 
Kwangju,  501-091 
Republic  of  Korea 

11.  Library 

P.O.box  77,  Gongneung-Dong 
Dobong-Gu,  Seoul,  132-240 
Republic  of  Korea 

12.  Library 

Army  Headquaters,  Bunam-Ri,  Duma-Myun 
Nonsan-Gun,  Chimgnam-Do,  320-919 
Republic  of  Korea 


